Astronomy & Astrophysics manuscript no. daSilvaetal2012 aph 


©ESO 2013 


March 2,2013 





Accurate and homogeneous abundance patterns in solar-type 

stars of the solar neighbourhood: a chemo-chronological analysis 

★ 

R. da Silva^ G.F. Porto de Mello^, A.C. Milone^ L. da Silva^ L.S. Ribeiro\ and H.J. Rocha-Pinto^ 



' INPE, Divisao de Astroflsica, Av. dos Astronautas, 1758 Sao Jose dos Campos, 12201-970 Brazil 
^ UFRJ, Observatorio do Valongo, Ladeira do Pedro Antonio, 43 Rio de Janeiro, 20080-090 Brazil 
' Observatorio Nacional, Rua Gal. Jose Cristino 77, Sao Cristovao, Rio de Janeiro, 20921-400 Brazil 

Received / accepted 

ABSTRACT 

Aims. We report the derivation of abundances of C, Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, Zr, Ba, Ce, Nd, and 
Sm in a sample of 25 solar-type stars of the solar neighbourhood, correlating the abundances with the stellar ages, kinematics, and 
orbital parameters. 

Metliods. The spectroscopic analysis, based on data of high resolution and high signal-to-noise ratio, was differential to the Sun and 
applied to atomic line equivalent widths supplemented by the spectral synthesis of C and C2 features. We also performed a statistical 
study by using the method of tree clustering analysis, searching for groups of stars sharing similar elemental abundance patterns. 
We derived the stellar parameters from various criteria, with average eiTors of 30 K, 0.13 dex, and 0.05 dex, respectively, for Ted, 
log^, and [Fe/H]. The average error of the [X/Fe] abundance ratios is 0.06 dex. Ages were derived from theoretical HR diagrams and 
membership of the stars in known kinematical moving groups. 

Results. We identified four stellar groups: one having, on average, over-solar abundances (([X/H]> = +0.26 dex), another with under- 
solar abundances (([X/H]> = -0.24 dex), and two with intermediate values (([X/H]) = -0.06 and +0.06 dex) but with distinct chemical 
patterns. Stars sharing solar metallicity, age, and Galactic orbit possibly have non-solar abundance ratios, a possible effect either of 
chemical heterogeneity in their natal clouds or migration. A trend of [Cu/Fe] with [Ba/Fe] seems to exist, in agreement with previous 
claims in the literature, and maybe also of [Sm/Fe] with [Ba/Fe]. No such congelation involving C, Na, Mn, and Zn is observed. The 
[X/Fe] ratios of various elements show significant correlations with age. [Mg/Fe], [Sc/Fe], and [Ti/Fe] increase with age. [Mn/Fe] 
and [Cu/Fe] display a more complex behaviour, first increasing towards younger stars up to the solar age, and then decreasing, a 
result we inteipret as possibly related to time-varying yields of SN la and the weak s-process in massive stars. The steepest negative 
age relation is due to [Ba/Fe], but only for stars younger than the Sun, and a similar though less significant behaviour is seen for Zr, 
Ce, and Nd. [Sr/Fe] and [Y/Fe] show a linearly increasing trend towards younger stars. The [Cu/Ba] and [Sm/Ba] therefore decrease 
for younger stars. We found that [Ba/Mg], [Ba/Zn], and [Sr,Y,Ba/Sm] increase but only for stars younger than the Sun, whereas the 
[Sr/Mg], [Y/Mg], [Sr/Zn], and [Y/Zn] ratios increase linearly towards younger stars over the whole age range. 

Key words, stars: solar-type - stars: fundamental parameters - stars: abundances 
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. 1. Introduction for numerous chemical elements, a successful model must har- 

rSl ' ™ „ , . , • , , , • , , • , „ ^ , monise stellar evolution inputs such as the initial mass func- 

^ . The Galactic chemical and dynamical hi story can be well framed ^-^^^ star-formation rate, and mass loss processes, connect- 

- : bLA|£nes oL|^rage laws (see e.g. [ Edvardsson et al. | l993 t ■ ^^^^^ -^^ properties of Galactic components, and 

.| McWilliam| | l997|) , namely: the age-metalhcity relation (which ^ j .^^ ^^^j .t^uctures, differing 

. in principle is accessible by more than one chemical element), fo^ stellar nucleosynthetic yields and their sensitiv- 

5_] . the [element/element] abundance ratios, the stellar metalhcity .^j^^ differing metallicities. 

rrt , frequency distribution, the Galactic metallicity gradient, and the ^ . , , . j- / ^ , • ■• , ■ • • 

. . . ^ „ ^- u- . 1. J • .II.- u- During the evolution of the Galactic disc, nucleosynthesis in 

star formation history, besides their mutual relationships as a . . ^ > • > ■ • , 

„ . J, J . rr,. r /-^ 1 u • successive generations of stars occurs together with dynamical 

function of space and time. The success of the Galactic chemi- . „ ^ , ^ , 

, ji-.u j ju.u- . J interactions with the interstellar gas. The states of the Galaxy 

cal evolution models IS to be ludged by their ability to reproduce . . , ^. , • ,1 111 

these constraints (see lMlen & Porto dTMdtolll^ . past periods of its evolution are still preserved in the abun- 

T , \, — r — r '■ — ^, ^ dance distributions of solar-type stars, which constitute an ideal 

In recent years there has been growing recognition that, even , . n,, , 

^, , , , ■ c^ AC A . 1 population to study the chemical evolution. These stars have an 

though such average laws are meaningful and fundamental, there ,. . , .c.^i t,, 

, ., ?, J 1 • 1 • 1 /-^ 1 age dispersion comparable to the age of the Galaxy. They are 

may be considerable underlying complexity in the real Galaxy . , ^ . , • , • > 

^, I , .1 . .1 • . J T . similar to the Sun in many physical parameters, allowing the ap- 

that has gone at least partly unappreciated. In the present age ,. . j-..^,. . . . 

c , . • J . {. J ■ u J plication of a differential analysis and the consequent minimisa- 
of very large spectroscopic databases and precise abundances \ ^ , .,, . _c , ■ , , , 
tion of theoretical shortcomings of atmospheric models and sys- 

Send offprint requests to: R. da Silva, tematic errors. In addition, their chemical composition does not 

e-mail: dasilvr2(agmail . com change in consequence of the mixing processes in their surfaces, 

* Based on observations collected at the Cerro Tololo Inter-American which means that the present abundance of a given element is the 

Observatory, Chile. same as in the time of their formation (an exception are the abun- 
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dances of Li, Be, and B, but these elements are not considered 
in this work). Therefore, the chemical abundance of solar-type 
stars, combined with kinematical, orbital, and evolutionary pa- 
rameters (mass and age), provide a powerful tool to investigate 
the chemical and dynamical evolution of the Galaxy. 

Over the last decade, several works have analys ed the com- 
position of disc dwarf sta r s of spectral types F and G ( Chen et al] 
2000^, 'Reddv et al.l l2003t iBensbv et alJ l2005t IChen et al.1 120081 



Neves et al. . 2009^ ). The metallicity ranges have become wider, 
the number of stars and elements studied has become larger, 
stars with and without detected planets have been compared, and 
the chemical distinction between thin and thick stellar popula- 
tion has been refined. Here we specifically ask, when regarded 
in as high a level of detail as possible with present techniques, to 
what extent the relative abundances of chemical elements can be 
traced as a function of age, the nature of underlying nucleosyn- 
thetic processes, and whether these properties can be statistically 
grouped for the nearby solar-type stars, defining "snapshots" rel- 
evant to the chemo-chronological evolution of the Galaxy. 

In this work we present a multi-elemental spectroscopic 
analysis of a sample of 25 solar- type stars in the solar neigh- 
bourhood, all members of the thin disc stellar population (ex- 
cepting one star in the transition between thin and thick discs). 
We have performed the determination of atmospheric parameters 
(effective temperature, metallicity, surface gravity, and micro- 
turbulence velocity), mass, age, kinematical and orbital param- 
eters, and elemental abundances based on equivalent widths or 
spectral synthesis. Three different criteria were used to pin down 
the stellar effective temperatures, and they showed excellent in- 
ternal consistency. We have also performed a statistical study of 
our a bundance results us ing the method of tree clustering anal- 
ysis ( lEveritt et al.l 1200 ll) . through which we looked for stellar 
groups that share similar abundances in [X/H], where X repre- 
sents one given element. Four groups were identified and then 
analysed in terms of their relations with [Fe/H], age, [Ba/Fe], 
kinematics, and Galactic orbits. 

Despite the small range in metallicity, our sample stars cover 
a broad range in age, and possible trends with age were traced. 
The relation between [X/Fe] and [Ba/Fe] for a few elements has 
also been considered given the Ba-rich nature of some of our 
stars. In particular, we investigated previous correlations with Na 
and Cu sugg ested by ICastro et al. (1999). Finally, considering 
the results of Rocha-P into et al.l (2006) that, on average, metal- 
poor and old stars tend to have larger |/?,n-/?0| (where is the 
mean orbital distance from the Galactic centre), we have looked 
for any relations involving the stellar groups of the clustering 
analysis and the kinematic and orbital parameters of the sample. 

Though limited in size, our sample was carefully built up to 
undergo an homogeneous and detailed analysis, based on spectra 
with high resolution and high signal-to-noise ratio, in order to 
achieve a precision as high as possible in our determinations. 
Particular care was exercised to derive the stellar atmospheric 
parameters from different and independent criteria, in an attempt 
to limit the abundance uncertainties as much as possible. 

In Sect. |2] we describe the observations and the reduction 
process. In Sect. |3]we present the methods used to derive the at- 
mospheric parameters and the chemical abundances. The stellar 
evolution, kinematics, and orbits are presented in Sect. 2] The 
tree clustering method is described in Sect. |5] and all the results 
are discussed in Sect. |6] Finally, we present our conclusions in 
SectH 



2. Observations and data reduction 

The sample stars wer e selected fr om the B right Star 
(lHoffleit& Jaschekll982h and Hipparcos (lESAll997h catalogues 
according to the following conditions: 

i) Solar neighbourhood stars in a distance < 40 pc; 

ii) Stars brighter than V = 6.5 and with declination < +20 °; 
Stars with eff'ective temperature and metallicity distributed 



Hi) 



iv) 



over about 5500 < reff< 6100 K and -0.3 < [Fe/H] 
< 0.3 dex, respectively, which represents the colour index 
range 0.52 < {B - V) < 0.78 from the {B - V) calibration 
described in Sect. l3.4l (Eq.[3Ti: and 

Stars with no information of duplicity (capable of sig- 
nificantly aff'ecting the spectroscopy) available in the 
astrometric and sp ectros co pic binary catalo gues of 



Hoffleit & Jascheg (l9S^ IWarren & Hoffleil 



Batten etal.l (Il989l) . and iDuquennov & Mavo3 (Il99lh : 



1987), 



the possibility of duplicity was a fterward revised in the 
Washington Double Star Catalogue (Maso n et ai]|2001h and 
in the survey of Raghavan et al. (2010), and no close-in 
companions that could affect our analysis were found. 

Based on these criteria, 99 F, G, and K dwarfs and subgiants 
were selected, out of which 25 stars were observed and analysed 
in this work. Our sample, shown in Table [T] contains F and G 
dwarf and subgiants stars from the thin disc stellar population, 
excepting the star HD 50806, which is probably in the transition 
between thin and thick discs. 

The observations were carried out at the Cerro Tololo Inter- 
American Observatory (CTIO, Chile) in two different runs: /) 15 
stars were observed in March 25-26, 1994 using the Cassegrain 
echelle spectrograph mounted on the 4 m telescope, with the 
red camera, 140 yum slit, and Tek CCD detector of 1024 x 1024 
pixels (24 X 24 yum pixel size), and with a gain of 1 e"/ADU; 
the spectra have resolution R ~ 29 000 and cover the wave- 
length range 4370-6870 A divided into 46 orders; and /;) 10 
stars were observed in November 8-15, 1997 using the bench- 
mounted echelle spectrograph and a 750 mm folded Schmidt 
camera attached to the 1.5 m telescope; the same CCD was used; 
the spectra have resolution R ~ A6 000 and cover the wavelength 
range 4550-6520 A divided into 37 orders. 

The two subsamples, although observed in different condi- 
tions, were both selected based on the same criteria and will be 
treated as an homogeneous single sample. The spectra collected 
in the first run have, on average, signal-to-noise ratio (S/N = 395 
± 60) slightly larger than those of the second one (S/N - 320 + 
90) and, despite having smaller resolution, may provide smaller 
uncertainties in some parameters estimated here. Nevertheless, 
all spectra have S/N > 200 in the blue region, which warranties 
spectral line profiles good enough to the equivalent width mea- 
surements. Any differences in the error estimates are discussed 
throughout the paper whenever needed. 

Two spectra of the sunlight reflected by Ganymede were also 
observed, one in each run. The S/N ratios were estimated using 
continuum windo ws in the spectra sel ected by inspection of the 
solar flux atlas of iKurucz et al.l ([1984^ (hereafter th e Solar Flux 
Atlas) and the solar line identifications catalogue of lMoore et al.l 
(1 19661) (hereafter the Solar Lines Catalogue). The mean val- 
ues ( S/N ) measured in the wavelength range 4500-5000 A 
are listed in Table [T] For larger wavelengths the S/N ratios 
are even higher, approaching twice that for /14500 in the range 
6000-6500 A. 
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Table 1. The 25 sample stars. The mean values of S/N ratios 
measured in continuum windows around AA 4500-5000 A are 
listed for the stars and for the two sunlight spectra (Ganymede). 



Object 


V 


Sp. type 


<S/N) 


First run (March 26-25, 1994): 


Oanymede 


^ 1 n 

J.IU 


LJZ V 






S 94. 


V 




HD 43834 


5.08 


G7 V 


370 


HD 84117 


4.93 


F8V 


410 


HD 102365 


4.89 


G2V 


370 


HD 112164 


5.89 


F9V 


400 


HD 114613 


4.85 


G3V 


400 


HD 115383 


5.19 


GOV 


320 


HD 115617 


4.74 


G7V 


370 


HD 117176 


4.97 


G5V 


410 


HD 128620 


-0.01 


G2V 


450 


HD 141004 


4.42 


GOV 


520 


HD 146233 


5.49 


G2V 


420 


HD 147513 


5.37 


G5V 


290 


HD 160691 


5.12 


G3 IV-V 


400 


HD 188376 


4.70 


G5 IV 


310 


Second mn (November 8-15, 1997): 


Ganymede 


5.10 


G2V 


410 


HD 1835 


6.39 


G3V 


230 


HD 26491 


6.37 


Gl V 


310 


HD 33021 


6.15 


GlIV 


240 


HD 39587 


4.39 


GOV 


510 


HD 50806 


6.05 


G5V 


300 


HD 53705 


5.56 


GOV 


370 


HD 177565 


6.15 


G6V 


290 


HD 181321 


6.48 


G2V 


220 


HD 189567 


6.07 


G2V 


260 


HD 196761 


6.36 


G8V 


400 



The spectra were reduced using routines for order 

identification and extraction, background subtraction (including 
bias and scattered light), flat-field correction, wavelength cali- 
bration, radial-velocity shift correction, and flux normalisation. 
The wavelength calibration was performed onto the stellar spec- 
tra themselves using lines selected by inspection of the Solar 
Flux Atlas and the Solar Lines Catalogue. The normalisation of 
the continuum is a very delicate and relevant step in the analy- 
sis procedure, since the accuracy of the equivalent width mea- 
surements is very sensitive to a faulty determination of the con- 
tinuum level. Therefore, continuum windows free from telluric 
or photospheric lines were carefully selected also based on the 
Solar Flux Atlas and the Solar Lines Catalogue. 



3. Spectroscopic analysis 

A differential spectroscopic analysis relative to the Sun was per- 
formed to determine the atmospheric parameters and the chem- 
ical abundance of several elements in our sample. The analysis 
was based on the equivalent widths of atomic lines measured in 
the spectra, and on the spectral synthesis of carbon atomic and 
molecular lines. The two groups of stars (15 observed in the first 
and 10 in the second run) were treated in comparison with the 
Ganymede spectrum of their respective run. 



' Image Reduction and Analysis Facility, distributed by the National 
Optical Astronomy Observatories (NOAO), USA. 



3.1. Equivalent widths and atomic line parameters 

Atomic lines of the elements Na, Mg, Si, Ca, Sc, Ti, V, Cr, Mn, 
Fe, Co, Ni, Cu, Zn, Sr, Y, Zr, Ba, Ce, Nd, and Sm were selected 
throughout the spectral range for equivalent width (EW) mea- 
surements. The lines were chosen based on the Solar Flux Atlas 
and the Solar Lines Catalogue, selecting only those for which 
the profiles were sufficiently clean from blends in order to pro- 
vide reliable measurements. For both solar and stellar spectra, 
the EW values of more than 7500 lines of these elements (about 
300 lines per star) were measured by hand by Gaussian function 
fit using IRAF routines. 

Strong line profiles are better described by Voigt functions 
than by Gaussian functions. We show in Fig. [1] (left panel) a 
comparison of the EWs measured in this work in the Ganymede 
spectrum of the second observation run (Gany 2) by Gaussian 
function fit to those measured in the Solar Flu x Atlas (R > 
500000 and S/N -3000) by Voigt function fit ( iMevlan et all 
Il993h . A similar diagram was obtained using the Ganymede 
spectrum of the first run (Gany 1) and the following relations 
represent the linear least square regressions fitted to both dia- 
grams: 

^'W'Solm- Flux Atlas = (0.3 + 0.6) + (1.065 ± 0.012) EWcanyl (1) 
^'W'Solm- Flux Atlas = (0.8 ± 0.5) + (1.075 ± 0.010) £WGany2 (2) 

where EW is given in mA. The standard deviations and the 
cross-correlation coefficients are, respectively, cr = 3.4 mA and 
r = 0.991 for Eq. [T] and cr = 2.7 mA and r = 0.994 for Eq. |2l 
Therefore, to reduce possible systematic uncertainties and pro- 
vide direct comparison with other works, all our EWs were 
transformed to a common system using the regression coeffi- 
cients of these equations (the constant terms have no statistical 
significance within 2cr). The regressions were derived in order 
to have a direct transformation to the Solar Flux Atlas system. 

We also did the converse (using -EWcanymede vs. -EWsolai- Flux Atlas 

diagrams) for comparison, and the resulting regressions are com- 
parable with those of Eq.[T]and|2] within Icr. 

The wavelength and lower-level excitation potential (x) of 
the atomic lines used were taken from the Solar Lines Catalogue. 
The oscillator strengths (gf) were computed using a solar model 
atmosphere applied to the EWs of Ganymede (converted using 
Eq.[T]and|2]i i n order to pr ovide the standard solar abundances of 
Greves se & Noelsl (Il993h . The adopted solar abundances are of 
course inconsequential in a differential analysis. 

The solar and stellar model atmospheres were computed 
with a code kindly supplied by Dr Monique Spite (Meudon 
Obse rvatory, Paris) that i nterpo lates the model-atmosphere grid 
from lEdvardsson et al.l (Il993|). We use d an updated version 
of the original code from Spitel (1 19671) . The fundamental at- 
mospheric parameters (effective temperature T^ff, metallicity 
[Fe/H], surface gravity \ogg, and micro-turbulence velocity ^) 
and the population ratio of helium and hydrogen atoms («He/nH) 
are taken as input. For the Sun, the adopted parameters are Tetf = 
5777 K, log g= 4.44,^= 1.3 kms"', nHe/«H = 0.1, and log eg = 
7.50 (the solar Fe abundance). 

The spectral lines used and their parameters are listed in 
Table [8] in which the EWs are the raw values (before the con- 
version). We do not list the EWs of the other stars but they are 
available upon request. 

The atomic lines of the elements Mg, Sc, V, Mn, Co, and 
Cu have important hyperfine struct ure (HFS). The ir gf values, 
listed in Table |9] were taken from ISteff^enl (Il98l and also re- 
vised according to the EWs of Ganymede and the standard solar 
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11 11 1 11 11 1 11 11 1 11 11 1 11 11 1 11 11 1 11 11 1 11 11 1 11 11 1 , 11 y/l , 

^'''solar Flux Altas = 3 + bEW^^^^^ 'i/ / 




- a = 0.8 ± 
b = 1.075 
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Fig. 1. Left panel: comparison of the EWs measured i n this work in the Ganymede spectrum (second observation run) by Gaussian 
function fit to those measured bv lMevlan et al.l (1 19931) in the Solar Flux Atlas by Voigt function fit; filled circles represent 145 lines 
of several elements; the linear regression (solid line) is the same expressed by Eq. |2] Right panel: comparison of the EWs listed in 
Table[8]of all lines in the Ganymede spectra in common to both runs after the conversions according to Eq.[T]and|2] the coefficients 
of the linear regression (solid line) are shown. 



abundances of iGrevesse &Noelsl(ll993h . as done for the other 
elements. For the elements Zn, Sr, Y, Zr, Ce, Ba, and Nd, for 
which theoretical HFS exist, either their effects are negligible or 
the spectral lines used are too weak to depend on the HFS as- 
sumpt ion. For the elements without HFS data listed in Ste ffenI 
(Il985h . values of neighbouring multiplets were adopted. The 
only exception was the Mg line at /15785.285, for which no HSF 
information was available. Its gf value, listed in Table [8] was 
obtained in the same way as those lines without HFS. This is 
not a strong line (ZsW < 60 mA) so that the eiT or induced in the 
Mg ab undance determination is not important. Idel Peloso et alJ 
(l2005h have recently shown that for Mn and Co the differences 
in their abundances computed using different values of HFS are 
not greater than 0. 10 dex. 

The HFS of Ba, and also its isotopic splitting, are of some 
importance only for the line at /16496.9, and can be neglected 
for the Unes at ^5853.7 and /i6141.7 (see Korotin et al. 20lil). 
However, for the 25 stars of our sample we have found a good 
agreement among the abundances yielded by the three lines, with 
a mean standard deviation of 0.07 dex. Moreover, a test per- 
formed using only /15853.7 and /16141.7 indicated that the global 
behaviour and trends found in the abundance diagrams, as well 
as all our conclusions involving the Ba results would not change 
if only these two lines were used. 

Concerning the fact that we used two spectrographs with 
different spectral resolutions, we performed a test in which we 
degraded the spectrum of Ganymede and of the metal-rich star 
HD 1835, both observed in the second observation run, to match 
the resolution of the first run. New values of EWs were then ob- 
tained, and no systematic differences were found when they are 
compared with the original measurements. Moreover, a compar- 
ison of the EWs of Ganymede listed in Table[8]and converted ac- 
cording to Eq.[T]and|2](see Fig.[T] light panel) demonstrates that 
the equivalent widths of the two observation runs were properly 
transformed to a common system, reinforcing our assumption of 
an homogeneous analysis. 



3.2. Derivation of atmosplieric parameters 

In order to determine the fundamental atmospheric parameters 
we developed a code that iteratively calculates these parameters 
for a given star based on initial input values. The so called exci- 
tation effective temperature (T^^'^) was calculated through the ex- 
citation equilibrium of neutral iron by removing any dependence 
in a [Fe 1/H] vs. x diagram. The micro-turbulence velocity was 
obtained by removing the dependence of [Fe 1/H] on EW, and 
the ionisation surface gravity (loggion) was computed through 
the ionisation equilibrium between Eel and Fell. Finally, the 
metallicity was yielded by the EW of Fe I lines. 

The temperature used in our abundance analyses was the ex- 
citation effective temperature, which is a better representation of 
the temperature stratification of the line forming layers. In order 
to compute the stellar luminosity with more reliability by also 
accounting for any consequence due to small LTE departures, 
we also considered two other temperature indicators, which are 
described in Sect. 13.41 

Table|2]lists the spectroscopic atmospheric parameters of the 
program stars. An estimate of their uncertainties was performed 
based on the analysis of HD 146233 and HD 26491, which are 
representative stars in our sample (from the first and second runs, 
respectively) with regard to their parameters and quality of the 
spectroscopic data. For both stars we have similar errors and they 
were obtained as follows: 

/) The uncertainty in T^^^ is related to the standard error of 
the angular coefficient of the linear regression fitted to the 
[Fe 1/H] vs. X diagram. The temperature is changed until this 
coefficient is of the same order of its error The difference 
between the best value and the previous one from the last 
iteration provides the uncertainty o-{T^^^) - 30 K; 

ii) The uncertainty in metallicity is the standard deviation of 
the abundance yielded by individual Eel lines, which is 
cr([Fe/H]) = 0.05 dex; 

///) To estimate the uncertainty in log^ion, its value is changed 
until the difference between the averaged abundance yielded 
by Fe 1 and Fe 11 lines is of the order of their internal er- 
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Table 2. Atmospheric parameters from our spectroscopic analy- 
sis. For comparison, the evolutionary surface gravity computed 
as described in Sect. l4.1l (Eq.l7Ti is also shown. 



Star logg.on logftvoi [kmVi] [Fe/H] 

0.00 
0.21 
0.22 
0.09 
0.20 
0.00 
0.11 
0.02 
0.22 
0.06 
0.28 
0.32 
0.15 
0.23 
0.00 
0.04 
0.23 
0.03 
0.05 
0.04 
0.27 
0.08 
0.06 
0.00 
0.27 
0.32 



rors (~0.05 dex), which led to an uncertainty cr(loggion) = 
0.13 dex; 

iv) The uncertainty in the micro-turbulence velocity is estimated 
regarding the [Fel/H] vs. EW diagram. The ^ value is mod- 
ified until the angular coefficient of the regression is of the 
same order of its error An uncertainty cr(^) = 0.04 kms"' 
was found for this parameter 

The spectroscopic atmospheric parameters were used to 
compute the model atmospheres, which in turn are required 
in the abundance determination. I n our analysis, we used the 
model-atmosphere grid derived by lEdvardsson et al.l (Il993b for 
stars with eff'ective temperatures from 5250 to 6000 K, sur- 
face gravity from 2.5 to 5.0 dex, and metallicity from -2.3 to 
-1-0.3 dex (with small extrapolations when needed). These are 
ID, plane-parallel, constant flux, line-blanketed, and LTE mod- 
els computed over 45 layers. 

The model atmospheres are, essentially, subject to errors 
in the atmospheric parameters, in the LTE simplifications, and 
in the thermal homogeneity assumption. However, the effects 
of non-LTE and thermal inhomogeneities are hopefully mi- 
nor for the elements and the stellar types studied here, being 
more important for low metallicity and low surface gravity stars 
(lEdvardsson et al.lll993l:lAsplundi2005h . Possible errors induced 
by such simplified assumptions are dominated by other sources 
of uncertainties. In addition, in a differential analysis, the errors 
in the stellar atmospheric structure are of second order 

We also investigated what would be the effects on the de- 
rived abundances if another model-atmosphere grid were used. 
We compared Edvardsson and Kurucz models and, using the 
same equivalent widths, gf values, and atmospheric parameters 
obtained from a solar spectrum, we found that the differences in 



Table 3. A comparison between the estimated abundance errors 
CTest (see Sect. 13.31 and I3.5l l and the dispersions of these abun- 
dances around the mean cr^isp for elements with N = 5 or more 
lines measured. For each element of each observation run the 
larger value was adopted to represent the errors cr([X/Fe]). 



HD 146233 HD 26491 







O" disp 


N 


0"esl 


O" disp 


N 


c 


0.07 


_ 


2 


0.07 


_ 


2 


Na 


0.03 


_ 


2 


0.06 


_ 


2 


Mg 


0.03 


_ 


4 


0.06 


_ 


4 


Si 


0.05 


0.06 


11 


0.06 


0.03 


17 


Ca 


0.03 


0.04 


6 


0.05 


0.05 


13 


Sc 


0.06 


0.03 


6 


0.09 


0.03 


13 


Ti 


0.07 


0.04 


24 


0.10 


0.04 


38 


V 


0.07 


0.07 


8 


0.11 


0.04 


11 


Cr 


0.06 


0.04 


14 


0.07 


0.04 


29 


Mn 


0.04 


0.05 


8 


0.06 


0.04 


11 


Co 


0.07 


0.05 


9 


0.11 


0.04 


12 


Ni 


0.04 


0.03 


23 


0.07 


0.04 


26 


Cu 


0.04 




3 


0.07 




3 


Zn 


0.05 




1 


0.06 




1 


Sr 


0.06 




1 


0.07 




1 


Y 


0.07 


0.05 


5 


0.09 


0.05 


5 


Zr 


0.09 




3 


0.08 




3 


Ba 


0.06 




3 


0.08 




3 


Ce 


0.07 


0.07 


5 


0.12 




3 


Nd 


0.16 




2 


0.12 




2 


Sm 


0.09 




1 


0.15 




1 



abundance [X/Fe] for most of the elements are of the order of 
0.03 dex or smaller, achieving a maximum of 0.07 dex. We note 
that these values are for the Sun and represent the differences 
between the two sets of model atmospheres. The effects of these 
differences when computing the stellar abundances relative to 
the Sun are minimised in a differential approach. 

3.3. Abundance determination and tineir uncertainties 

The abundance of the elements studied were determined using an 
adapted version of a code also supplied by Dr Monique Spite. 
The code takes into account the solar gf values and the stellar 
model atmospheres (computed using the atmospheric parame- 
ters of each star) to calculate the abundances that fits the equiv- 
alent widths measured in the spectra (transformed according to 
the procedure described in Sect. 13.1) . The results of this abun- 
dance determination are presented and discussed in Sect. |6] 

The main sources of uncertainties in the abundance determi- 
nation come from the errors in the EWs (the most important), the 
gf values, the atmospheric parameters, and the adopted model 
atmospheres (these two latter are discussed in Sect. 13.2) . 

The uncertainties in the EWs, were estimated as follows: by 
plotting the diagram /JWho 26491 vs. -EWcanymede and computing 
the standard deviation of the linear regression, we obtained cr - 
2.9 mA. The solar EWs were measured in the surrogate spectrum 
of the Sun collected under the same circumstances as for the 
program stars. Therefore, we assumed that cr is a quadratic sum 
of the errors in EW of both objects and that they are similar to 
each other. Thus, for the star HD 26491 the value of criEW) is 
cr/ V2 = 2.1 mA. Similarly, for the star HD 146233 we obtained 
cr{EW) = 1.7 mA. These values were adopted to represent the 
uncertainties in EW of the two observation runs. Because the 
solar gf values were computed to reproduce solar the equivalent 
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Table 4. Colour indices and photometric effective temperatur es (given in K). The values of ( B - V) and (Bj - Vt) were taken from 
the Hipparcos Catalogue . The Ha effective temperatures from lLvra & Porto de Mellol(l2005l) are also listed. For the star HD 128620 
(a Cen A), the photometric determination of Tf.f{ was not performed (see discussion in Sect. 13.4b . 
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widths, the errors in EW contribute twice to the total uncertainty, 
with approximately the same magnitude. 

Each one in turn, EW, T^^'^, [Fe/H], ^, and log^ion are 
changed by Icr in the sense of increasing the abundance ratios 
and new abundances are computed for each element. The dif- 
ferences between new and previous abundance values provide 
the errors induced by each parameter and a quadratic sum of 
these errors yields the total uncertainty in the elemental abun- 
dance ratios. The estimated errors (aest) are listed in Table[3]for 
both HD 146233 and HD 26491, and they are compared to the 
dispersions around the mean (cr^isp) for elements with at least 
five lines measured in the spectra. For these elements, the larger 
values were adopted to represent the errors cr([X/Fe]) in each 
observation run. Otherwise, o-g^i was adopted. 

3.4. Photometric and Ha effective temperatures 

The effective temperature of the sample stars were also obtained 
using some photometric calibrations, providing the photomet- 
ric effective tempera t ure (T ^g"')- These calibrations, derived by 
iPorto de Mello et al.l (|20 l ife r the {B - V), {Bj - Vj), {b - y), 
and p colour indices, are given by the following equations; 

^eff(B-y) = 7747 - 3016 [B -V)[\- 0.15 [Fe/H]} (3) 
^eCx-Vx) = 7551 - 2406 (Bt - Vt) {l - 0.2 [Fe/H]} (4) 
T^T(t-y) = 8481 - 6516 {b -y)[\- 0.09 [Fe/H]} (5) 



^phot 



= 11654 yS- 2.349 



(6) 



for Teff given in K. The standard deviations of these calibrations 
are cr = 65, 64, 55, and 70 K, respectively. 

The {B - V) and {Bj - Vt) colour indices of our stars were 
taken from the Hipparcos Catalogue, and {b — y) and fi from the 
literature (see Table |4]i, when availa ble. For the star H D 33021, 
the p values adopted are only from iPerrv et all (Il987b because 
these authors made 41 measurements of this index. For the star 
HD 50806, only one reference for the [} index was found, and 
the effective temperature from this index strongly disagrees with 
that obtained from the other colours and we thus discarded it. 

Table |4] lists the colour indices used and the photometric 
effective temperature (T^^) obtained, which is a mean of the 
temperatures computed using the four calibrations, weighted by 
their variances. The referenc es for {b — y) and (i are als o hsted. 
The v a lues of [ b - y ) from Gronbech & Olsen (Il976h . lOlsenI 
(ll983h . lTwarog (1980), and lSchuster & Nissenl(ll988l) were con- 
verted to the Olsen (19931) system according to equations pro- 
vided by the latter author 

The star HD 128620 is the primary component {V - -0.01) 
of the a Cen triple system. The Teff determination for very bright 
stars using photometric colours is normally considered risky due 
to systematic effects that may affect the results (such as non- 
linearity and detector dead time) and also, in the case of this 
system, due to a possible contamination by the companion. For 
this reason, we preferred do not include this star in our T^s esti- 
mates based on the photometric indices. Nevertheless, our spec- 
troscopic determination for a Cen A, - 5857 + 30 K, is in 
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Fig. 2. Spectral synthesis of the C2 molecular band regions (/15128 and /15165) and of the C atomic lines (/i5052.2 and /i5380.3) for 
Ganymede of the second observation run (top panels) and for one metal-rich star of the first run, HD 160691 (bottom panels). The 
solid lines represent three models with diff'erent values of [C/Fe]: -0.05 (red), 0.0 (green), and +0.05 dex (blue). The differences 
between observed (dots) and computed spectra (o-c) are also shown in the bottom of each panel. 



go od agreement with the phot ometric determination performed 
bv lPorto de Mello et al.l (l2008h . T^s = 5794 + 34 K. 

Concerning the uncertainties in on the one hand, the 

internal error of the weighted mean, computed using the stan- 
dard deviations in the four photometric calibrations, is 31 K. 
On the other hand, the mean value of the standard deviations 
around (i.e., the dispersion of the four values of temper- 
ature around the weighted mean) is 41 K. Therefore, in this 
work we adopted cr{T^^) = 40 K as the internal uncertainty in 
our photometric effective temperatures. This uncertainty intrin- 
sically takes into account the errors in the colour indices them- 
selves. 

The stellar effective temperature can also be estimated by 
modelling the wing profile of th e Ha line, whi ch is very sensi- 
tive to changes in this parameter iLvra & Porto de Mello (2005|) 
applied this method to solar neighbourhood stars, and values of 
for our sample were used as a third T^s indicator, with an 
uncertainty criT"^^) = 50 K. 

3.5. Carbon abundance from spectral synthesis 

The carbon abundance was derived using the spectral synthesis 
method applied to molecular lines of electronic-vibrational band 
heads of the C2 Swan System at /15128 and /15165, and also to 
C atomic lines at /15052.2 and /15380.3. To reproduce the atomic 
and molecular absorption lines in the observed spectra of the 
sample stars, the MOOG spectral synthesis cod^E developed by 
Chris Sneden (University of Texas, USA), was used. The syn- 
thetic spectra were computed in steps of 0.02 A, also taking into 
account the continuum opacity contribution in ranges of 0.5 A. 



^ http://www.as.utexas.edu/~chris/moog.html 



The Unsold approximation multiplied by 6.3 was adopted in the 
calculations of the line damping parameters. 

The model atmospheres are the same used in the spectro- 
scopic analysis. They also include the micro-turbulence velocity 
and the elemental abundances of each star, both assumed to be 
constant in all layers. For any element X for which no abundance 
was determined in this work, we adopted the metallicity of the 
respective star to set the [X/H] ratio. 

The atomic and molecular line parameters used to compute 
the synthetic spectra are: the central wavelength, the gf values, 
the lower-level excitation potential, and the constant of dissocia- 
tion energy Dq (only for molecular features). Atomic and molec- 
ular data were taken, respectively, from the Vienna Atom ic Line 
Database - VALD (Kupkaetal. 1999, 2000; Piskun ov et alj 
I1995I: iRvabchikova et al.Hl997h and from|K uruczl (II992I) . In ad- 
dition to C2, the spectral regions studied also include MgH 
molecular features that may contribute to the continuum for- 
mation. The oscillator strengths of C2 and MgH lines were re- 
vised according to the norm alisation of the Honl-London factors 
(IWhiting & Nichollsl[T97l . 

To account for the spectral line broadening, the synthetic 
spectra were computed by means of the convolution with three 
input parameters: /) the spectroscopic instrumental broadening; 
//) the limb darkening of the stellar disc; and ///) a composite of 
velocity fields, such as rotation velocity and macro-turbulence 
broadening, named Vbioad- The instrumental broadening was esti- 
mated by means of the FWHM of thorium lines present in Th-Ar 
spectra observed at the CTIO. The linear limb-darkening coeffi- 
cient (on average u ~ 0.7 for all the sample stars) was individ- 
ua lly estimated by interpolatio n of and log gio„ in Table 1 
of lDfaz-Cordovez et al.l (fT995h . As a first estimate of Vbroad, the 
projected rotation velocity (vsin /) of the stars was used, which 
was computed based on the profile of four isolated Fe I lines 
(/15852.2, /15855.1, /15856.1, and /15859.6) in the spectra. Small 



7 



R. da Silva et al.: Chemo-chronological analysis of solar-type stars 



^ 200 



CUV 

I 

v% _200 



I 

t_ -200 



I 



-200 



nean ± (7 = -10 ± 42 




-i — I — I — I — \ — I — I — I — \ — I — I — I — I — I — I — \ — |- 

mean ± cr = 25 ± 46 



"A" 
i 




Tmean 
^ eff 



[K] 



Fig. 3. Comparisons of the three effective temperature indica- 
tors: the excitation vs. photometric (top panel), excitation vs. Ha 
(middle panel), and photometric vs. Ha (bottom panel) tempera- 
tures. The symbols represent the stellar groups defined in Sect.|5] 
(see Fig.|6]l. 

corrections in Vbioad were applied when needed according to an 
eye-trained inspection of the synthetic spectra. The final values 
are listed in Table |5] where they can be compared to the stellar 
age and the chromospheric activity level. 

Figure |2] shows two examples of synthetic spectra of the C2 
molecular band regions around /15128 and /15165, and of the 
C atomic lines at /15052.2 and /15380.3 for the sunlight spec- 
tram reflected by Ganymede (second observation ran) and for 
the metal-rich star HD 160691. The spectral synthesis was first 
applied to the Ganymede spectra of both runs, then the gf values 
of some atomic and molecular lines were revised when needed, 
and finally the synthesis was applied to the other stars, treated 
according to their observation rans. For each case, the best fit 
was obtained through the minimisation of the nns between ob- 
served and synthetic spectra. 

In order to estimate the uncertainties in the [C/Fe] abundance 
determination, we performed a spectral synthesis of the most 
prominent molecular band used (/15165) adopting model atmo- 
spheres perturbed by the errors estimated for the atmospheric pa- 
rameters. This procedure resulted in: +0.03 dex due to the error 
in r^^.'=; ±0.01 dex due to the error in [Fe/H]; +0.02 dex due to 
the error in ^; and ±0.03 dex due to the error in log ^io,,. The un- 
certainties related to errors in Vbioad (~ 1 -0 km s"' or smaller) and 
in the limb darkening coefficient are negligible. The quadratic 
sum of the individual contributions (also including a global er- 
ror of 0.05 dex estimated based on the nns minimisation of the 
solar spectrum) yields a total uncertainty cr([C/Fe]) - 0.07 dex. 



4. Evolutionary, Kinematic, and orbital parameters 

4.1. Mass and age determination 

The evolutionary parameters mass and age were obtained by 
interpolation in the Yonsei-Yale (Y^) evolutionary tracks and 



isochrones (lYi et al.llIOOll iKim etalJl2002h drawn on the HR 
diagram and computed for difi'erent values of metallicity. 

The luminosity used in the diagrams were calculated with 
parallaxes taken from the new reduction of the Hipparc os data 
(va n Leeuwenll200"7h . bolometric corrections (BC) from iFlowerl 
(19961), and an absolute bolometric magnitude for the Sun 
M® , - 4.75 for M® - 4.82. We estimated, for these nearby 
stars with precise parallaxes, a mean uncertainty of 0.01 dex in 
log(L/L0). 

The effective temperature is a weighted mean (T™^") of the 
excitation, photometric, and Ha temperatures. The weights are 
given by for cr(r«^'=) = 30 K, cr(rP^^= 40 K, and o-{T^p 
- 50 K, obtained as described in Sect. 13.21 and 13.41 The uncer- 
tainty of the weighted mean is 22 K, calculated using these three 
values of cr. On the other hand, the standard deviations of the 
three values of temperature around the weighted mean is 29 K. 
Therefore, we conclude that our estimates of effective temper- 
ature based on the three indicators agree with each other very 
well (see the comparison in Fig.[3]l and that the mean value has 
a mean internal error cr{T'^'^^) = 30 K. 

iPorto de Mello et all (|2008|) determined the effective temper- 
ature of Q-Cen A (HD 128620) and B also using the excitation, 
photometric, and Ho- approaches. They found a good agreement 
for a Cen A, a solar temperature star. However, for a Cen B (reft 
~5200 K), the excitation effective temperature is about 100- 
150 K higher than the photometric and the Ha counterparts, 
which the authors attributed to non-LTE effects. Although the 
agreement for the coolest and hottest stars in our sample is not 
that good, especially in the comparison of excitation and Ha 
temperatures. Fig. |3] does not show any systematic difference 
among the three indicators, and the differences are nonetheless 
within 2cr for all the Tgft range (only HD 1 15383 has 7'™'^ larger 
than by 3cr). This confirms that ID LTE model atmospheres 
may adequately represent solar-type stars, at least in a differen- 
tial analysis relative to the Sun. 

We remind that the temperature used in our abundance anal- 
yses (based on the equivalent widths or synthesis of spectral 
features) was the excitation effective temperature, which better 
characterises the temperature radial profile in the stellar pho- 
tosphere and the formation of absorption lines in the emergent 
spectrum. On the other hand, to better represent the luminosity 
of a star and to account for any possible effect due to small de- 
viations from LTE, we adopted the weighted mean of the three 
temperature indicators. 

The stars were grouped according to their values of metallic- 
ity (12 groups from [Fe/H] = -0.32 to +0.32 dex) and then their 
masses and ages were computed using evolutionary tracks and 
isochrones for each stellar group. The difference in metallicity 
between each star and its respective HR diagram is not greater 
than 0.02 dex. A few examples for some metallicities are shown 
in Fig.|4] To reproduce the Sun's position in the diagrams, adopt- 
ing T gff - 5777 K and age = 4.53 Gy dGuenther & Demarquel 
Il997h . the evolutionary tracks and isochrones were displaced in 
log(reff) and log(L/L0) by 0.001628 (-22 K in Teft) and 0.011, 
respectively. These values are, at any rate, of the same order or 
smaller than the uncertainties on these parameters. 

As an independent check, we calculated the evolutionary sur- 
face gravity using the values of mass and effective temperature 
obtained, which we called gsmi, using the following equation: 



^ y^mean 

log ^ = log — + 4 log + 0.4(Mboi - Ml,) 



(7) 
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Fig. 4. Evolutionary tracks (magenta solid lines) and isochrones (blue dots and dashed lines) from Yonsei-Yale, calculated for 
different values of metallicity, showing how to derive the stellar masses and ages. The ages in Gyr are indicated alongside the 
evolutionary tracks, for which the masses are indicated in Mq. The Sun's position for T^g - 5777 K is also shown in the solar 
metallicity panel (top-left corner). 



where Mboi is the absolute bolometric magnitude for the stars. 
The values of log ^evoi are listed in Table|2]together with the ion- 
isation surface gravity. They are in very good agreement, having 
a dispersion of only 0.09 dex, smaller than the uncertainty of 
0.13 dex estimated for loggion- 



4.2. Galactic velocities, distance, and eccentricity 

The kinematic properties of our sample were investigated 
by computing the Galactic velocity components J/lsr, ^lsr, 
and Wlsr (see Fig. |5]l with respect to the Local Standard 
of Rest (LSR). We d e velope d a code that uses equations of 
iJohnson & SoderblomI (Il987h . parallaxes and proper motions 
both from the new reduction of the H i pparcos data, and ra dial 
velocities from iHolmberg etaLl (l2007h . iTorres et al.l (|2006|) for 



HD 114613, and Santos et al 



2004 for HD 160691. For the 
Sun, the adopted values of Ulsr, Vlsr, and Wlsr are 10.0, 5.3, 
and 7.2 kms respectively dDehnen & Binnevlll998l) . 

The mean orbital distance from the Galactic centre (Rm) and 
the orbital eccentricity (e) were also considered in our analysis 
(see Fig.|5j, where e = (R^- Rp)/(R^+ Rp) and R^ = (R^+ Rp)/2 
were computed using the perigalactic (Rp) and the apogalac- 
tic {R„) orbital distanc es from the Geneva-Copenhagen survey 
(IHolmberg et al.ll2009l) . For the Sun, the adopted values are e - 
0.06 and R^^S kpc. 



5. Tree clustering analysis 

We looked for statistically significant abundance groups in our 
sample using a hierarchical clustering analysis. To avoid missing 
abundance values, the analysis uses only those elements having 
abundances measured for all stars and was applied to the [X/H] 
abundance space. 

We us ed the complete lin kage method for the hierarchical 
clustering (lEverittetal.ll200lh and euclidean distances as mea- 
sures of dissimilarities in this abundance space. A hierarchical 
clustering algorithm works by joining similar objects in a hier- 
archical structure. Initially, each object is assigned to its own 
cluster The algorithm proceeds iteratively, joining the two most 
similar clusters in each pass until there is just a single cluster 
The resulting hierarchy of clusters for our data is shown in Fig.|6] 
(upper panel) as a dendrogram. In this plot, the most similar ob- 
jects are linked together in the bottom forming clusters, which 
are then iteratively linked together in pairs by similarity. The 
vertical axis in a dendrogram measures the dissimilarity between 
each individual or cluster Since we used euclidean distances in 
the [X/H] abundance space, the units of this axis is dex, although 
it measures the total dissimilarity in this abundance space and 
not in a single variable. 

Clusters can be defined by specifying a reasonable total dis- 
similarity value for pruning the dendrogram. There is no un- 
ambiguous or optimal way for defining this pruning value, es- 
pecially because the clusters found depend on the clustering 
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Fig. 5. Galactic velocity components (left and middle panels) and orbital eccentricity as a function of the mean orbital distance from 
the Galactic centre (right panel) for the sample stars and the Sun. The symbols represent the stellar groups defined in Sect. |5] 



method and cluster shapes. Since our sample is quite small, we 
arbitrarily decided to prune our dendrogram at the total dissimi- 
larity of 0.9 dex, which is shown in Fig.|6]by the dashed red line. 
This pruning value was chosen in order to have four more or less 
equally populated clusters. Considering the size of our sample, 
less than four clusters would simply limit our discussion to poor 
against rich stars, while a larger number of clusters would make 
such an analysis meaningless. 

The average ([X/H]) behaviour of the four clusters for each 
element considered in this analysis is shown in the middle panel 
of Fig.|6l Two of these clusters have over-solar abundances, with 
averages +0.26 and +0.06 dex for all the elements, whereas 
the two others have under-solar abundance values, with aver- 
ages -0.06 and -0.24 dex. We can observe in the figure that the 
elemental abundance patterns of the metal-poor and metal-rich 
groups are distinct between each other In particular, it seems 
to exit a chemical distinction even between the two intermedi- 
ate groups: the under-solar intermediate group has an abundance 
pattern that roughly follows the element by element pattern of 
the metal-poor group, whereas the pattern of the over-solar in- 
termediate group resembles the scaled-solar mixture. 

The clustering analysis we have presented was tentatively 
based on biological ideas of evolving species, in the sense that 
the material the stars came from is continuously changing. For 
this reason we made use of [X/H] abundances ratios instead of 
[X/Fe]. We implicitly need the time evolution that [X/H] has, 
because we want a time hierarchy in the output groups. A sim- 
ilar analysis in the [X/Fe] space can still show groups, but the 
hierarchical relation between these groups in a dendrogram will 
not necessarily show evolutionary trends, because this variable 
is only indirectly linked with time. Notwithstanding, we have 
checked this, but the output groups show no meaningful interpre- 
tation in terms of chemical evolution or abundance ratio groups. 
The outcome could be different if the sample were larger, but 
this needs to be verified with another sample, what is beyond the 
scope of this paper Nevertheless, we included in Fig. |6] (bottom 
panel) the average ([X/Fe]) behaviour of the stars clustered in 
the [X/H] parameter space. The small variation of ([X/Fe]) with 
respect to the solar values reinforces our point above: that for 
this specific small sample, the [X/Fe] parameter space is dynam- 
ically very narrow and does not favour a cluster analysis. 



6. Results and discussion 

Table|5]lists the evolutionary (mass and age), kinematic {Ulsr, 
Vlsr, and Wlsr velocities), and orbital (mean orbital distance 
from the Galactic centre and orbital eccentricity) parameters 
computed for the program stars. They are grouped following the 
tree clustering analysis performed in Sect. |5] 

The uncertainties in the mass and age determination may 
vary widely depending on the stellar position in the HR diagram. 
We made an estimate of these errors for a few representative stars 
in our sample (cool and hot dwarfs and subgiants). We took into 
account the errors estimated for log(L/L0), log(reif), and also 
[Fe/H] considering that the evolutionary tracks and isochrones 
are metallicity dependent. We found that the uncertainties in 
mass stand between 0.02 and 0.08 Mq, whereas those in age vary 
from about 0.5 Gyr (or smaller) for evolved stars up to about 
2.5 Gyr for cool main-sequence stars. 

For HD 1835, HD 39587, HD 147513, and HD 181321 we 
determined an approximative value for their masses and an 
upper limit for their ages given their position in the HR 
diagram (close to the Zero Age Main Sequence). Indeed, 
these are very young stars: one of them, HD 1835, is likely 
a member of the Hyades s t ar clu ster (~600 Myr) accord- 



ing to iLopez-Santiago et alJ ilOld): two o thers, HD 39587 



(Soderblom & Mayor" 1993; Fuhrmann"2004) and HD 147513 
tSoderblom & Mavoriil993: .Montes et al.n2001ft b elong to the 
Ursa Major moving group of ~300 Myr (see also ICastro et al.l 
IT999I) : and HD 181321 is a m ember o f the Castor moving group 
(~200 Myr) according to Montes et al. (2001). For these four 
stars, we adopted the ages of their respective moving group in 
our study. 

The stars HD 1 12164 and HD 160691, indicated by asterisks 
(*) in Table |5] are located in a region of the HR diagram where 
successive evolutionary tracks and isochrones are superposed 
(see example in Fig.|4|. Therefore, their mass and age determi- 
nation may yield larger uncertainties: < 0.12 M© and < 0.7 Gyr 
for HD 1 12164, and < 0.04 Mq and < 1.0 Gyr for HD 160691. 

The C/lsr, V'lsr, and Wlsr velocities have a typical un- 
certainty of 0.3 kms"' or smaller An exception is the star 
HD 188376, for which the large errors (2.5, 1.1, and 1.8 kms ^ 
respectively) are due to a large uncertainty in its parallax. The 
level of activity in the chromosphere of the stars, which is related 
to their age, was also investigated. The t able Usts the flux in the 
centre of the Ha line (F' ), computed bv lLvra & Porto de Mellol 
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(I2005h and used as a chromospheric activity indicator (the larger 
the value of F^^, the higher the level of chromospheric activity). 
The uncertainty for this parameter is 0.5 x 10^ erg cm"^ s"' . 

We note that a few stars in our sample have at least one 
planetary companion detected. They are HD 102365, HDl 15617, 
HDl 17176, HD147513, and HD160691 (see The Extrasolar 
Planets Encyclopaedia: http://exoplanet.eu). Comparisons of 
properties of stars with and without planets are frequently pub- 
lished. In the analysis performed in this paper, however, no pe- 



cuUar information distinguishing the two populations has been 
found. 



6.1. Elemental abundances 

Table [TOl lists the chemical abundances relative to iron [X/Fe] 
obtained for the elements studied. For some elements of some 
stars, the abundance determination was not possible due to the 
poor quality or the weakness of their spectral lines (empty fields 
in the table). The carbon abundance ratios [C/Fe] are shown in 
Tabled 

We also computed the mean abundance ratios ([X/Fe]) of 
the following groups of elements: ;) two groups of light metals: 
(Mg, Si) and (Ca, Sc, Ti); //) two groups of the iron peak: (V, 
Cr, Co, Ni) and (Mn, Cu); Hi) light elements from the s-process: 
(Sr, Y, Zr), to which we refer as Is; and iv) heavy elements from 
the s-process: (Ba, Ce, Nd), referred to as hs (see Table |7). They 
were grouped either because they have possibly the same nucle- 
osynthetic origin or because they share a similar behaviour in 
the diagrams. The abundance ratio between heavy and light el- 
ements from the s-process, [hs/ls] = [hs/Fe] - [Is/Fe], was also 
calculated. 

Figure [7] shows diagrams with the abundance ratios of 
the program stars for individual elements and nucleosynthetic 
groups. The uncertainties are listed in Tables |3] and |7] As for 
the individual elements, the estimated errors are compared to 
the dispersions around the mean for groups having at least two 
elements. For each group of each observation run (represented 
by the stai's HD 146233 and HD 26491), the larger values were 
adopted to be the uncertainties in the grouped abundance ratios. 
In Fig.|2l the stars are represented by different symbols accord- 
ing to the tree clustering analysis performed in Sect. |5] 

The star HD 1835 is enriched in Ca, Sr, and Ba. The mean 
value of all s-process elements also suggests an over-solar abun- 
dance. Sm, the only r-process element analysed here, shows an 
under-solar abundance of -0.3 dex, but with a large error. As al- 
ready mentioned, this is a very young star, a probable member of 
the Hyades star cluster of age ~600 Myr, which is in agreement 
with its high level of chromospheric activity indicated by F'^^ in 
Tableia 

Two other very young stars are HD 39587 and HD 147513, 
both members of the kinematic Ursa Major group. They are 
clearly overabundant in the s-process elements, especially Ba, 
and unde rabundant in C, which is in agr eeme nt with the 
results of iPorto de Mello & da Silval (Il997ah and ICastro et al" 



. These stars were proposed bv lPorto de Mello & da Silva 



I to be barium stars, originated in a phenomenon in which 
the more massive component of a binary system evolves as a 
thermally pulsing asymptotic giant branch (TP-AGB) star and 
the material produced in the He-burning envelope, enriched in s- 
process elements, is dredged-up to the surface and then accreted 
by its companion by wind mass transfer The initially more mas- 
sive star is now a white dwarf whereas the companion has be- 
come the primary barium star At present, HD 39587 is a single- 
lined spectroscopic and astrometric binary, with a low mass 
companion of 0.15 Mq (Konig et al. 2002), and HD 147513 has 
a common proper motion co mpanion, a DA2 whi te dwarf, at an 
angular separation of 345" dHolberg et al. 2002|). Th e barium- 
star scenario was not supported by Castro et al.l ( 19991) . who pro- 
posed that the two stars simply have usual Ba abundances for 
their age and that probably all the Ursa Major group members 
are Ba-enriched, either due to a primordial origin or because they 
are young (see discussion in Sect. l6.4l i. 
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Table 5. Ev olutionary, kinematic, and orbit al parameters, separating the stars according to the classification of Sect.|5] The Ha fluxes 
(^Hff)' from lLvra & Porto de Mellol (12005V are given in units of 10^ ergcm"^ s The mean distance from the Galactic centre (Rm) 
and the orbital eccentricity (e) are from Holmberg et al..(,2009i) . The asterisk (*) indicates stars adjacent to superposed evolutionary 
tracks and isochrones on the HR diagram. 
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HD 181321 and HD 188376 are two other Ba-rich stars. 
HD 188376 is the most evolved and massive star analysed here, 
clearly in the evolutionary stage of a subgiant. The other star, 
HD 181321, is the youngest and has the highest level of chro- 
mospheric activity in our sample. Indeed, our determination for 
Vbroad is 12.5 kms"', indicating a fast-rotating star. It has so- 
lar atmospheric parameters, excepting a high value of micro- 
turbulence velocity = 2.3 kms '). The kinematic and orbital 
parameters are also very close to the solar values. In other words, 
this star has, on the one hand, about the same effective tempera- 
ture, metallicity, surface gravity, mass. Galactic orbit, and space 
velocities as the Sun. On the other hand, it is very young and 
significantly enriched in Ba, strengthening the relation between 
Ba abundance and age (see Sect. 16.4b . 

The high microturbulence velocity of the star HD 181321 is 
probably prompted by the strengthened convection and turbu- 
lence in its upper photosphere, which is subjected to large non- 
thermal energy influxes from the chromosphere. The UV radia- 
tion excess from the chromosphere of an active star can scape 
to the photosphere and cause departures from LTE due to an 
ionisation imbalance. The induced overionisation is commonly 
manifested by differences either in excitation and photometric 
effective tem peratures, or in ionisation and evolutionary sii rface 
gravities (see lPorto de Mello et al.ll2008HRibas et al. Il2010l) . For 



HD 181321, our determination of and T^^ are in very good 
agreement with each other Therefore, only the difference in sur- 
face gravity and the large value of microturbulence velocity are 
possible signs that an overionisation is taking place in the pho- 
tosphere of this active star A full non-LTE analysis and a photo- 
spheric and chromospheric modelling would probably settle the 
issue, but this goes beyond the scope of this paper 

It is also very worthwhile to investigate the chemical abun- 
dances in stars that share similar values of age, metallicity, and 
Galactic orbit {R^ and e). These subgroup of stars are sup- 
posed to share the same physical conditions of the Galaxy at 
the time and galactocentric position of their birth. An example 
of this includes the stars HD 43834, HD 841 17, HD 141004, and 
HD 146233, which also have ages, metallicities, and Galactic or- 
bits close to the solar values (they were all classified as inter- 
mediate abundance stars in the tree clustering analysis, one in 
slight underabundance and the three others in slight overabun- 
dance with respect to the Sun). In spite of this, only HD 43834 
and HD 141004 show solar abundances, within the uncertain- 
ties, for all (or almost all) elements. HD84117 is deficient in 
Mn and enriched in Na and in eleme nts of the s-process (Sr, Y, 
Zr, and Ba). HD 146233. proposed bv lPorto de Mello & da Silval 
(Il997bh as the closest solar twin ever known at that time, is ac- 
tually (as also proposed by these authors) enriched in some el- 
ements of the s-process (Sr and Ba) and possibly enriched in 
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Fig. 7. Abundance ratios [X/Fe] for individual elements (top of each panel) and ([X/Fe]) for the nucleosynthetic groups (bottom of 
each panel). The uncertainties adopted are those listed in Table s|3] and [7] The symbols represent the stellar classification defined in 
Sect. |5](see Fig. ill. 
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Fig. 7. continued. 



Sc, V, and Sm. Thus, possibly, the investigation of a larger sam- 
ple of stars at a similar level of detail as done here could reveal 
that non-solar abundance ratios are present even for stars shar- 
ing the same place and time of birth. Whether this reflects intrin- 
sic heterogeneities in their natal interstellar clouds (in principle 
a reasonable hypothesis since the elements reflecting non-solar 
ratios are related to different nucleosynthetic processes operat- 
ing in different timescales) or else is evidence for considerable 
radial migration in the Galaxy is a question we plan to address 
in a subsequent work involving a larger sample. 

Still concerning the relations involving the stellar groups 
from the clustering analysis and the kinematic and orbital pa- 
rameters of our sample, we can see in Fig. |5]that the group of 



metal-poor and old stars seems to have larger velocities in the 
direction of the Galactic centre (IL^lsrI > 40 kms"') and larger 
eccentricities (e > 0.2) than the other stars. The star HD 50806 
appears to have a singular position in this figure (in particular, 
it has the most eccentric orbit among the sample stars), which 
is probably related to its membership in the transition popula- 
tion of thin-thick disc stars. The limitation of our sample does 
not allow to verify the results of iRocha-Pinto et aH (l2006l) that 
metal-poor and old stars show more orbital radial spread in the 
Galaxy. 
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Table 6. Carbon abundance ratios derived from the spectral syn- 
thesis of the two atomic lines (/15052.2 and /15380.3) and the two 
C2 band heads (/15128 and /15165) studied. Column 6 shows the 
final values adopted. 



Table 7. Groups of chemical elements and uncertainties in the 
mean abundance ratios ([X/Fe]). The estimated errors are com- 
pared to the dispersions around the mean. For each group, the 
larger value of each observation run was adopted. 



Star 


[C/Fe] 


[C/Fe] 


[C/Fe] 


[C/Fe] 


rC/Fel 


/i5052 


/{5380 


.15128 


/!5165 


HD 1835 


0.07 


0.05 


0.02 


0.04 


0.05 


HD 20807 


0.06 


0.05 


-0.11 


0.05 


0.01 


HD 26491 


0.01 


0.01 


-0.06 


-0.04 


-0.02 


HD 33021 


-0.05 


-0.04 


-0.14 


-0.19 


-0.11 


HD 39587 


-0.07 


-0.07 


-0.12 


-0.12 


-0.10 


HD 43834 


-0.04 


0.01 


0.01 


0.02 


0.00 


HD 50806 


0.21 


0.13 


0.08 


0.08 


0.13 


HD 53705 


-0.07 


-0.04 


-0.14 


-0.07 


-0.08 


HD 841 17 


0.09 


0.07 


-0.06 


-0.06 


0.01 


HD 102365 


0.07 


0.12 


-0.08 


0.02 


0.03 


HD 112164 


-0.01 


-0.01 


0.00 


0.02 


0.00 


HD 114613 


-0.01 


-0.02 


0.03 


0.00 


0.00 


HD 115383 


-0.12 


-0.09 


-0.05 


-0.05 


-0.08 


HD 115617 


-0.11 


-0.04 


0.00 


0.04 


-0.03 


HD 117176 


-0.06 


-0.04 


-0.06 


-0.04 


-0.05 


HD 128620 


0.01 


0.04 


0.03 


0.07 


0.04 


HD 141004 


-0.03 


-0.03 


0.00 


0.07 


0.00 


HD 146233 


-0.04 


-0.05 


-0.07 


0.01 


-0.04 


HD 147513 


-0.15 


-0.11 


-0.18 


-0.10 


-0.14 


HD 160691 


0.01 


-0.01 


0.02 


0.06 


0.02 


HD 177565 


0.05 


0.10 


0.02 


0.03 


0.05 


HD 181321 




-0.10 


-0.15 


-0.15 


-0.13 


HD 188376 


-0.16 


-0.06 


-0.02 


-0.05 


-0.07 


HD 189567 


0.00 


-0.05 


-0.15 


-0.10 


-0.08 


HD 196761 




-0.13 


-0.20 


-0.13 


-0.15 



6.2. Abundance trends as a function of[Fe/l-l] 

Through the analysis of Fig. |8]we investigate possible trends in 
the abundance ratios as a function of the stellar metallicity. These 
trends are more clearly identified if the elements are grouped to- 
gether, either based on their nucleosynthetic origin or because 
they share similar trends in the diagrams. For this reason, we 
show in the bottom panels of this figure the mean abundance 
{[X/Fe]) of a few groups of elements as a function of the metal- 
licity (the same groups plotted in Fig. |7] and listed in Table |7]i. 
The stars are also identified according to the tree clustering anal- 
ysis. We have fitted linear regressions on the diagrams in three 
ranges of metallicity: for stars poorer than the Sun, for stars of 
solar metallicity or richer, and for all the sample stars. We have 
then computed the cross-correlation coefficients in these three 
metallicity ranges and plotted the regressions of the more signif- 
icant trends (only if | r | > 0.5). 

The overall trend of our abundance ratios as a function 
of the stellar metallicity normally follows what has been sug- 
gested in the literature concerning the nucleosynthet ic origin of 
the elements and their abundance evolution in time ( Chen et al.! 
2000tlReddvetal.l l2003t iBensbv et al.ll2005t ICherTet al., .20081: 
Neves et al I l2009h . The light metals Ca, Sc, and Ti are predomi- 
nantly produced by Type II Supernovae (SN II) at the beginning 
of the enrichment history of the Galactic disc. On the other hand, 
iron and the iron-peak elements V, Cr, Co, and Ni are predomi- 
nantly synthesised by Type la Supernovae (SN la) in longer time 
scales. Therefore, it is expected that the abundance ratio of these 
light metals with respect to iron progressively decrease from 
metal-poor to metal-rich stars, whereas the abundance ratio of 
iron-peak elements with respect to iron, all produced at the same 
rate, remains constant and close to zero in the whole range of 
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Ba, Ce, Nd 


0.06 
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0.06 


0.07 



metallicity. Indeed, this is exactly what is observed in Fig.[8]for 
[Ca, Sc, Ti/Fe] and [V, Cr, Co, Ni/Fe], to within our stated abun- 
dance uncertainties. 

The light metals IVIg and Si may not only be produced 
by SN II considering that [Mg, Si/Fe] flattens out for metallic- 
ities higher than -0. 1 dex. The same beh aviour was found by 
IChen et all (l2000h and iNeves et all (l2009h . who suggested that 
SN la is possibly contributing. The star HD 50806 is clearly en- 
riched in Mg (and in other elements as well) according to Fig. [8] 
probably reflecting its membership to the thin-thick disc transi- 
tion. 

The situation of Mn, Cu, and Zn is somewhat more complex. 
The hypothesis of production in SNIa still stands, but this i s 
probably not the unique source. [Alien & Porto de Mellol (1201 ll) . 
in their study of s-process enriched stars, suggested that SN la is 
the main sourc e of production of ma nganese, in opposition to the 
conclusions o f iFeltzing eta"l] ( l2007h . who suggested that this el- 
ement is mainly produced by SN II. Our results in Fig. |8] which 
show an increasing trend of [Mn/Fe] as a function of [Fe/H] (see 
also the bottom panel of Fig. |6] in which there is a sequential 
crescent ordination of {[X/Fe]) from the metal -poor clustering 
group to the metal-rich one), seem to support the idea of an ex- 
tra nucleosynthetic source for the Mn yields. Such an increasing 
trend is also usually attributed to a metallicity dependence in 
the production of Mn in both SN la and SN II. In this work, Mn 
and Cu were plotted together through the mean abundance ratio 
[Mn, Cu/Fe]. Both these elements have abundances that increase 
with metallicity, though for Cu this trend is not as significant as 
for Mn, and seems to happen only for higher metallicities, being 
constant and close to zero for [Fe/H] < 0. Cu and Zn, although 
being adjacent elements in the periodic table, stand in the tran- 
sition between iron-peak and s-process elements, and their be- 
haviour is in sharp contrast. A decreasing trend in [Zn/Fe] vs. 
[Fe/H] i s seen for stars poorer than the Sun, in agreement with 
Fig. 1 of lAUen & Porto de Mellol (|20TTI) . 

The elements C and Na were not grouped together with other 
elements. C is synthesised in several different sites and behaves 
similarly to O and N, with [C/Fe] decreasing with increasing 
metallicity. This negative trend is mostly observed in the metal- 
poor regime ([Fe/H] < -0.3), hence not seen in Fig. [8] given 
the limited metallicity range of our program stars. Nevertheless, 
our results agree ver y well with the recent C abundance determi- 
nation performed by Ida Silva et al ] (120111) for solar-like dwarfs. 
Na is probably produced, among other processes, in the core of 
massive stars and ejected by SN II into the interstellar medium. 
Here we found that [Na/Fe] is constant and nearly close to zero 
in the range of metal-p oor stars, wi t h a po ssible increasing trend 
for higher metallicities. IChen et al.l(l2000l) suggested that maybe 
[Na/Fe] is close to zero for the whole metallicity range of disc 
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Fig. 8. Abundance ratios as a function of the stellar metallicity for individual elements (top panels) and for nucleosynthetic groups 
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the 95% confidence intervals (hashed area), and the cross-correlation coefficients are also shown for | r | > 0.5. The symbols follow 
the classification of Sect.|5](see Fig.|6|. 
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Fig. 9. Abundance ratios as a function of [Ba/Fe]. For the elements of the light s-process (Sr, Y, and Zr) the mean abundance ratios 
are plotted. The linear regressions (solid line), the 95% confidence intervals (hashed area), and the cross-correlation coefficients are 
also shown for | r | > 0.5. The symbols follow the classification of Sect.|5](see Fig.|6]l. 



stars. iNeves et"an (l2009h . however, found that [Na/Fe] is close 
to zero for thin disc stars for [Fe/H] between -0.2 and -1-0.2, but 
above solar for other values of metallicity. 

The [X/Fe] abundance ratios of elements of the s-process, 
mainly produced in TP-AGB of intermediate or low mass stars, 
and the r-process, produced in sites with high neutron density 
such as the final stages of massive stars (SN II, neutron stars), are 
supposed, respectively, to progressively increase and decrease 
from metal-poor stars to higher metallicities. These facts reflect 
the production of these elements in different time scales with re- 
spect to iron, the former products of long-lived AGE stars, the 
latter arising from short-lived massive stars which explode as su- 
pernovae. This is not clearly observed in the diagrams of Fig. [8] 
for Sm and the s-process elements because of the short metallic- 
ity range. An interesting result of the analysis of these diagrams 
are the properties involving some Ba-enriched stars, which we 
discuss in the next section. 

6.3. Abundance trends as a function of [Ba/Fe] 

A few stars in our sample are remarkably enriched in Ba by 
more than 3o-, especially HD 39587 and HD 147513. For this 
reason we also investigated the behaviour of the abundance ra- 
tios [X/Fe] or ([X/Fe]) as a function of [Ba/Fe] for some ele- 
ments or groups of elements showing some kind of relation with 
the production of barium (see Fig.|9]l. Once more, we computed 
the cross-correlations coefficients and we plotted in the figure 
the regressions of the most significant trends (| r | > 0.5). 

ICastro et al.l (1 1999b proposed the existence of an anticor- 
relation between the abundances of Cu and the s-process ele- 
ments. They found that [Cu/Fe] decreases with the increasing 
of [Ba/Fe], suggesting a relation between the destruction of Cu 
and the production of Ba (and other s-process elements) . In the 
recent analysis of Ba-enriched stars of lAllen & Porto de Melld 
( 1201 Ih . the authors have not supported this scenario, arguing that 
Cu seems to be little (or not at all) affected by the s-process. 



even though they have acknowledged that some Ba-rich stars do 
present anticorrelated abundances of Cu and the s-process el- 
ements. Our results point to a statistically significant decrease 
in th e abundances of Cu with increasing [Ba/Fe], in accordance 
with ICastro et al.l (Il999h . 

Two other iron-peak elements, Mn and Zn, are also shown in 
Fig.|9]and no clear correlation is observed. This may indicate that 
Mn and Cu do not share t he sam e nucleosynthetic origin. Indeed, 
lAllen & Porto de Mellol (1201 lb found that the synthesis of Cu re- 
ceives a larger contribution from n ot so massive stars than Zn, a 
result roughly in line with those of lCastro et al.l(ll999l) and ours. 
Such results point towards the necessity of both more extensive 
observations of the abundances of Cu and Zn, and more pro- 
tracted theoretical efforts, in order that a better understanding 
of the complex chemical history of these two elements may be 
achieved. 

ICastro et al.l (1 19991) also proposed an anticorrelation in the 
abundances of C and Na with respect to [Ba/Fe]. Our results do 
not seem to support this anticorrelation, though our most Ba-rich 
stars, the Ursa Major group members HD39587 and Hpi475I3 , 
are markedly C-deficient. iPorto de Mello & da Silval (Il997al) 
attributed the [C/Fe] deficiency of a barium star to the 
'^C(q', «)'*'0 reaction that occurred in the hot-bottom envelope 
of its companion during the TP-AGB phase. However, HD39587 
and HD 1475 13 are no longer regarded as true barium stars. 

Figure |9] shows an evident expected correlation for the light 
s-process elements (Sr, Y, and Zr). Correlations involving Nd, 
another heavy element of the s-process may also exist, but its 
abundance determination has larger uncertainties. An anticorre- 
lation between [Sm/Fe] and [Ba/Fe] also seems to exist. Sm is 
a good representative of the r-process elements, and despite the 
very large uncertainties in the determination of such elements, 
usually showing very few lines in the spectra of solar-type stars, 
an interpretation in which this anticorrelation is due to ever more 
efficient production of s-process elements in AGB stars, as com- 
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pared to the production of the r-process in SNII, seems war- 
ranted. 



6.4. Abundance trends as a function of age 

One of the motivations of the present paper is to explore the 
abundance ratios of elements due to different nucleosynthetic 
processes and the stellar ages, taking advantage of the reason- 
ably precise ages that can be attributed to our program stars. 
Figure [To] shows the diagrams [X/Fe] and ([X/Fe]) as a func- 
tion of the stellar age, and Fig. [TT] explores the relation be- 
tween the abundance ratios [X1/X2] of different elements and 
age. Similarly to Fig.|8] we have fitted linear regressions on the 
diagrams in three ranges of stellar age: for stars younger than the 
Sun (age < 4.53 Gyr), for stars older than the Sun, and for all 
the sample stars. We have then computed the cross-correlation 
coefficients in these three ranges of stellar age and plotted the 
regressions of the more significant trends (only if | r | > 0.5). In 
these figures there are positive, negative, or flat abundance trends 
in the three ranges of age. We notice, however, that the age of the 
Sun, used as a reference, was arbitrarily chosen. The exact value 
of the transition age when the abundance behaviour changes is 
not clear from these plots (it is a value between 4 and 6 Gyr). 

In spite of the long recognition (though not undisputed) of 
the so-called age-metallicity relation (see Fig.fTZb. the individual 
abundances in Fig. [TOland those in Fig. |8]do not share exactly 
the same behaviour, leading us to suggest in the following that 
the age-metallicity relation may be a multidimensional concept. 
For this reason, we regrouped the elements according to their 
abundance behaviour with age (see the bottom panels of Fig.fTOli. 

Carbon and sodium do not seem to present any important 
trends of [X/Fe] with age. The positive trends observed for Mg, 
Sc, and Ti (less clearly seen for Si and Ca) are simply the result 
of the Galactic chemical evolution. The production rate of these 
elements by SN 11 decrease with time since the formation of the 
Galactic disc (equivalently to increasing with age) compared to 
the increased production of Fe by the longer-lived SN la as we 
approach more recent epochs. Silicon perhaps shows no trend at 
all; Mg seems to have a more or less positive linear trend with 
increasing age; in their turn Ca, Sc, and Ti sport a more complex 
behaviour. The statistical significance of the behaviour of Ca is 
slight, and not much confidence should be placed in the appar- 
ent [Ca/Fe] decrease with time, followed by an increase towards 
more recent times. Taken at face value, this would appear to lend 
support to the suggestion that some fraction of the Ca synthesis 
might be due to SN la, in unison with their production of Fe. Sc 
and Ti seem to have a significant decrease with time with respect 
to Fe, but this decrease stops at a time close to the solar age and 
flattens thereafter towards present times. 

Among the Fe-peak elements, no important trend is seen 
in the [X/Fe] relation with age for V, Cr, Co, and Zn; only a 
dubious one for Ni in the interval of young stars. Yet, again, 
Cu and Mn suggest more underlying complexity. Even though 
the statistical significance of the linear regressions is slight, the 
abundance ratios to Fe of both these elements seem first to in- 
crease towards the present epoch, and then decrease (a behaviour 
that is reinforced when these two elements are plotted together 
through the mean abundance r atio [Mn, Cu/Fe]). Recalling that 
lAllen & Porto de Melld (1201 ll) have found that Mn is mostly 
due to SN la, our result could imply that the relative yield of Mn 
to Fe in SN la decreases with time (and consequently the over- 
all me tallicity). The situation for Cu is less straightforward, as 
usual. lAllen & Porto de Mellol ( 1201 ll) suggest that litde of the 
synthesis of Mn, Cu, and Zn is owed to the main s-process. 
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Fig. 12. Stellar metallicity as a function of age. The symbols fol- 
low the classification of Sect.|5](see Fig.|6]i. 



leaving the action of AGB stars an unlikely source of such a 
behaviour These same authors assert that the action of the so- 
called weak s-process, sited at the He-burning core of massive 
stars, has a non negligible contribution to the synthesis of Mn, 
Cu, and Zn. One possible explanation for the decrease of the 
abundance ratios of [Mn/Fe] and [Cu/Fe] towards more recent 
times is a decreasing yield of the weak s-process in their synthe- 
sis, as contrasted to the production of Fe by SN la. Clearly, an 
understanding of the detailed chemical evolution of elements, in 
both the dimensions of metallicity and age, of the Fe-peak and 
its transition with the heavier elements deserves closer scrutiny, 
both observationally and theoretically. 

We next turn to the [X/Fe] -age relation for t he s-process 
elements. These present particular in terest, since ICastro et al.l 



(Il999h suggested, also using data from lEdvardsson et al.l(ll993ir 
that [Ba/H] has a steeper decrease with age than [Fe/H], and 
therefore that [Ba/F e] increases towards modern times (see also 
iBensbv et al]l2007h . This fact might be interpreted, again, as a 
consequence of the larger yield of the s-process element synthe- 
sis by the long-lived AGB stars in relation to the not as-long- 
lived production of Fe by SN la. Do the other s-process elements 
show a behaviour similar to barium? Apparently this is so, but 
not in a straightforward way. The lighter s-process elements Sr 
and Y seem to have a linear trend of [X/Fe] with age (clearly 
seen for Sr), increasing towards present times, as expected. The 
[Zr/Fe], [Ce/Fe], and [Nd/Fe] ratio, however, appears flat in the 
old age regime, possibly increasing only for stars younger than 
the Sun. The [Ba/Fe] ratio behaves similarly but the increase to- 
wards younger ages is much sharper and more significant. At 
face value, these results point towards the evolution of the rel- 
ative yields of such s-process elements with time (and metal- 
licity) in AGB stars, apparently favouring the heavier species 
Ba and Nd over lighter ones. The simultaneous analysis of the 
[Cu/Fe] and [s-process/Fe] ratios, epitomised, for example, by 
the [Ba/Fe] relation with age, suggests an anticorrelation of Cu 
and Ba towards younger stars, as found bv ICastro et al.l (Il999l) . 
adding to the controversy surrounding the chemical evolution of 
copper. 

Finally, a positive trend with increasing age is observed 
for [Sm/Fe] in stars younger than the Sun, again reflecting the 
smaller number of SN 11 than SN la in the present in compari- 
son with the past (in this case the epoch of the Sun's formation). 
Significant positive trends are also observed for Ni and Cu in the 
regime of younger stars. 
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Fig. 10. Abundance ratios as a function of the stellar age for individual elements (top panels) and for nucleosynthetic groups (bottom 
panels). The vertical dashed line indicates the adopted solar age (4.53 Gyr). The linear regressions (solid line), the 95% confidence 
intervals (hashed area), and the cross-correlation coefficients are also shown for | r | > 0.5. The symbols follow the classification of 
Sect.|5](seeFig.|6|. 
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Fig. 11. Abundance ratios as a function of the stellar age. The vertical dashed line indicates the adopted solar age (4.53 Gyr). The 
linear regressions (solid line), the 95% confidence intervals (hashed area), and the cross-correlation coefficients are also shown for 
I r I > 0.5. The symbols follow the classification of Sect.|5](see Fig.|6]l. 



Concerning the four groups yielded by the tree clustering 
analysis and their relation with age, we can only state that, as ex- 
pected, the old stars in our sample are also metal-poor, whereas 
young stars tend to be metal-rich. Once more, the case of the 
star HD 50806 is evidenced. Although classified in the interme- 
diate group of stars with slightly over-solar abundances (([X/H]) 
= +0.06 dex), it is situated close to the group of metal-poor stars 
(([X/H]) — -0.24 dex) in Fig.[TOl which is maybe a consequence 
of its population membership. 

The bottom panels of Fig.[TO]shows the average of the [X/Fe] 
relations with age for selected groups of elements. Thus we see 
that grouping Sc and Ti reinforces the relation with age already 
shown by each element individually, and the same is seen for the 
grouping of Sr and Y. Similarly, averaging the [X/Fe] relations 
with age for Zr, Ba, and Nd produces a very steep increase to- 
wards present epochs, after a flat relation from the birth of the 
Galactic disc up to the solar age. 

These results prompted us to investigate the specifics of 
the [X1/X2] ratios with age, where Xi and X2 designate ele- 
ments other than Fe. In Fig.[TT]we explore the time evolution of 
some elements that displayed a particular clear [X/Fe] relation 
with age. Thus, it is apparent that the [Ba/Mg] and [Ba/Zn] in- 
crease steeply towards present epochs, for stars younger than the 
Sun, reinforcing their individual and opposite behaviour in the 
[X/Fe]-age diagrams. Also, both [Sr/Mg] and [Y/Mg] increase 
Unearly and significantly from the oldest to the youngest stars; 
the same is seen in the [Sr/Zn] and [Y/Zn] ratios. Investigating 
the [X1/X2] ratios of s-process to r-process elements, we found 
steep increases in the [Ba/Sm], [Sr/Sm], and [Y/Sm] ratios to- 
wards younger stars, but only for objects younger than the Sun. 
Significant age relations are, therefore, evidenced in the [X1/X2] 
ratios of diverse elements, representing a wide range of nucle- 
osynthetic processes and tentatively allowing the proposition 
that the age-metallicity relation is a more complex constraint to 
Galactic chemo-dynamical models than hitherto recognised. 



6.5. Abundance trends with condensation temperature 

Our determination of multi-elemental abundances also provides 
the study of possible trends in the abundance ratios [X/Fe] as a 
function of the condensation temperature (Tc) of each element. 
Values of 50% Tq (the temperature when 50% of an element is 
in the cond ensed phase) for a so lar-system compositi on gas were 
taken from lLoddersl(l2003l) and lLodders et al.l (120091) . 

For a few stars in our sample we have found some correla- 
tions of [X/Fe] with Tq (even after corrections due to Galactic 
chemical evolution effects were applied). However, no clear cor- 
relation seems to exist when comparing the slopes for refractory 
elements (those with Tc > 900 K) with several stellar parameters 
(r™™, [Fe/H], \ogg, ^, mas s, and age). The one involving the 
metallicity was proposed by iRamirez et al.l (|20I0|) in the sense 
that higher-metallicity stars present more negative slopes. Our 
results agree with their paper, but the number of metal-rich stars 
in our sample is too small to confirm their conclusions. 

7. Conclusions 

In this work we have performed a multi-elemental, differentially 
with respect to the Sun, spectroscopic analysis of a sample of 
25 solar-type stars in the solar neighbourhood. We have derived 
their atmospheric parameters (from various nearly independent 
criteria and with low internal errors), masses, ages, kinematical 
and orbital parameters, and elemental abundances (derived with 
very low internal uncertainties) based on equivalent widths or 
spectral synthesis. 

Despite small in size, our sample was carefully selected to 
undergo an homogeneous and detailed analysis based on spectra 
with high resolution and high signal-to-noise ratio. We have: (i) 
checked the effective temperatures based on the excitation equi- 
librium of neutral iron lines against those from photometric cal- 
ibrations and from the Ha wings profile, ( ii) checked the surface 
gravities computed through the ionisation equilibrium between 
Fe I and Fe II lines against those computed based on the evolu- 
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tionary parameters, ( Hi) derived masses and ages from evolution- 
ary tracks and isochrones computed considering the metallicity 
of each star, and (iv) appUed a differential spectroscopic analy- 
sis relative to the Sun, hence minimising the systematic errors 
and yielding a mean uncertainty of 0.06 dex in the abundance 
ratios. We thus expect that our determinations have achieved a 
high level of precision and accuracy. 

We have also applied a statistical study to our abundance re- 
sults using the method of tree clustering analysis, through which 
we looked for groupings of stars that share similar abundances 
in the [X/H] space. Although our sample has a limited range in 
metallicity, it covers a broad range in age. The detailed abun- 
dance pattern was then investigated through correlations with 
kinematics. Galactic orbits, and stellar ages. Our conclusions are 
thus summarised: 



1- Four groups were identified, two having over-solar abun- 
dances (with averages +0.26 and +0.06 dex on [X/H]), and 
two with under-solar abundance values (on average -0.06 
and -0.24 dex). Possible non-solar abundance ratios, even 
for stars which share the same age. Galactic orbit, and metal- 
licity as the Sun, are suggested. Whether these are due to 
heterogeneity in the stellar natal clouds, or by dynamical mi- 
gration within the Galactic disc, should b e investigated with 
larger samples. In particular, the results of lRocha-Pinto et aU 
(l2006l) that metal-poor and old stars show more orbital radial 
spread in the Galaxy could not be verified given the limita- 
tion of our sample; 

2- The presence of Ba-enriched stars in our sample prompted us 
to investigate in detail the relation of some elements with Ba. 
An anticorrelation between [Cu/Fe] and [Ba/Fe] was found, 
in line with similar claims in the literature. The [Sm/Fe] 
abundance ratios seem to be anticorrelated with [Ba/Fe], 
barely at the 95% confidence level. On the other hand, previ- 
ous suggestions of [C/Fe] and [Na/Fe] anticorrelations with 
[Ba/Fe] could not be confirmed. Even though the possible 
connected chemical evolution of Mn, Cu, and Zn has been 
recently discussed in the literature, no significant trend of 
[Mn/Fe] and [Zn/Fe] with [Ba/Fe] is suggested; 

3- The consideration of the [X/Fe] ratios with age revealed 
much differing behaviour of the elements, suggesting that 
the age-metallicity relation has more underlying complex- 
ity than commonly recognised. The [Mg/Fe], [Sc/Fe], and 
[Ti/Fe] decrease towards younger ages. The [Cu/Fe] and 
[Mn/Fe] ratios initially increase towards younger stars up to 
the solar age, and then decrease towards the youngest ob- 
jects, a result that may speak of differing yields in the SNIa 
production of these elements, related to metallicity and age, 
as well as a possible influence of the evolution with time 
of the weak s-process yields, operating off of massive stars. 
The steepest relation with age was found for the [Ba/Fe] ra- 
tio, but only for ages younger than the solar one, and a sim- 
ilar but less evident behaviour is seen for Zr, Ce, and Nd. 
Other heavy s-process elements, however, such as Sr and Y, 
show a linearly increasing [X/Fe] towards younger ages, par- 
ticularly clearer for Sr [Sm/Fe] significantly decreases for 
stars younger than the Sun. Thus, the [Cu/Ba] ratio clearly 
decreases towards younger stars, and the same is seen at a 
significant level for [Sm/Ba]; 

4- Considering the average of elements with similar behaviour 
with age considerably reinforces the aforementioned results, 
particularly for the [Sc,Ti/Fe], [Mn,Cu/Fe], [Sr,Y/Fe], and 
[Zr,Ba,Nd/Fe] relations; 



5- The consideration of element ratios not directly involving 
Fe shows some marked behaviour Particularly, the [Ba/Mg], 
[Ba/Zn], [Ba/Sm], [Sr/Sm], and [Y/Sm] steeply increase to- 
wards younger ages for stars younger than the Sun. Also, the 
[Sr/Mg] and [Y/Mg] ratios linearly increase towards younger 
ages, and the same is seen for [Sr/Zn] and [Y/Zn]; 

6- Possible correlations of the abundances, the condensation 
temperatures of the different elements, and the presence of 
exoplanets in our program stars was deeply investigated, but 
no significant correlation was found. 

The detailed consideration of precise element abundances 
derived from high-quality atmospheric parameters and spectro- 
scopic data, tied to masses, kinematics, and ages for solar-type 
stars, generally provides a wealth of interesting data, contribut- 
ing towards a broader understanding of the evolution of the 
Galaxy in its chemical and dynamical aspects. 
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Table 8. Atomic line parameters of the elements used in the analysis. Oscillator strengths (log gf) and raw EWs (before the conver- 
sion set out by Eq.[T]and |2]l, given in mA, of both the Ganymede spectra observed in the first (Gany 1) and second (Gany 2) runs 
are listed (except for C, for which the analysis is based on spectral synthesis). Lines with missing gf values represent the elements 
with hyperfine structure (Mg, Sc, V, Mn, Co, and Cu) and the detailed line splitting is shown in Table |9l 



^[A] 


Id. 


y 
[eV] 


Gany 1 


Gany 2 


A [A] 


Id. 


y 
[eV] 


Gany 1 


Gany 2 






t W 


log 8.f 


c W 




log gf 


c W 


log^f/ 


EW 


5052.167 


CI 


7.68 


-1.48 




-1.62 




4926.147 


Til 


0.82 






-2.17 


7.0 


5380 322 


CI 


7.68 


-1.78 




-1.84 




5022.871 


Til 


0.83 


-0.52 


113 


-0.35 


79.2 


6154.230 


Nal 


2.10 


-1.52 


42.3 


-1.52 


40.8 


5024.842 


Til 


0.82 


-0.66 


71.5 


-0.48 


73.8 


6160 753 


Nal 


2.10 


-1.23 


61.8 


-1.29 


59.9 


5071.472 


Til 


1.46 






-0.77 


31.3 


4571.102 


Ms I 


0.00 




107.8 




110.3 


5113.448 


Til 


1.44 


-0.86 


30.3 


-0.88 


27.6 


4730.038 


Ms I 


4.34 




78.2 




70.9 


5145.464 


Til 


1.46 


-0.67 


38.8 


-0.64 


37.6 


5711.095 




4.34 




119.2 




106.4 


5147.479 


Til 


0.00 






-1.98 


42.4 


5785 285 


Mai 


5.11 


-1.87 


50.5 


-1.82 


55.3 


5152 185 


Til 


0.02 






-2.03 


39.1 


5517 533 


Si I 


5.08 


-2.42 


15.7 


-2.51 


13.0 


5192.969 


Til 


0.02 


-1.10 


88.3 


-1.02 


84.4 


5621 607 


Si I 


5.08 






-2.61 


10.5 


5211.206 


Til 


0.84 


-2.10 


9.1 


-2.07 


9.2 


5665 563 


Si I 


4.92 


-1.98 


42.2 


-1.96 


41.8 


5219.700 


Til 


0.02 


-2.32 


27.9 


-2.23 


30.0 


5684 484 


Si I 


4.95 






-1.60 


62.3 


5295 780 


Til 


1.07 






-1.60 


13.2 


5690 433 


Si I 


4.93 


-1.80 


52.9 


-1.81 


50.6 


5426 236 


Til 


0.02 






-2.97 


7.8 


5701 108 


Si I 


4.93 


-1.90 


46.3 


-1.97 


41.0 


5471.197 


Til 


1.44 






-1.48 


9.6 


5708 405 


Si I 


4.95 






-1.35 


79.1 


5490 150 


Til 


1.46 






-1.00 


22.6 


5753 622 


Si I 


5.61 






-1.24 


50.8 


5648 567 


Til 


2.49 


-0.40 


11.5 


-0.39 


11.1 


5772.149 


Si I 


5.08 






-1.56 


57.7 


5679 937 


Til 


2.47 


-0.65 


6.9 


-0.63 


7.0 


5793 080 


Si I 


4.93 


-1.91 


46.0 


-1.94 


42.6 


5739.464 


Til 


2.25 


-0.67 


10.5 


-0.75 


8.6 


6125.021 


Si I 


5.61 






-1.50 


34.1 


5866 452 


Til 


1.07 






-0.82 


49.6 


6131.577 


Si I 


5.61 


-1.67 


25.8 


-1.65 


26.6 


6064 629 


Til 


1.05 






-1.88 


9.7 


6131 858 


Si I 


5.61 


-1.66 


26.6 


-1.64 


27.0 


6091.177 


Til 


2.27 






-0.44 


15.9 


6142.494 


Si I 


5.62 


-1.44 


37.5 


-1.45 


36.3 


6092 798 


Til 


1.89 


-1.31 


6.0 


-1.28 


6.1 


6145 020 


Si I 


5.61 


-1.40 


40.8 


-1.36 


41.8 


6098 694 


Til 


3.06 






-0.16 


6.1 


6243 823 


Si I 


5.61 


-1.22 


52.5 


-1.19 


52.7 


6126.224 


Til 


1.07 






-1.40 


23.5 


6244.476 


Si I 


5.61 


-1.26 


49.5 


-1.25 


48.6 


6258 104 


Til 


1.44 


-1.46 


54.6 


-0.43 


52.3 


5261.708 


Cal 


2.52 






-0.65 


99.4 


4568 345 


Till 


1.22 


-2.85 


33.6 


-2.85 


32.3 


5581 979 


Cal 


2.52 






-0.68 


97.2 


4583.415 


Till 


1.16 


-2.85 


36.2 


-2.84 


35.1 


5590 126 


Cal 


2.52 


-0.78 


96.3 


-0.73 


93.9 


4657 209 


Till 


1.24 






-2.31 


55.8 


5867 572 


Cal 


2.93 


-1.62 


25.3 


-1.59 


25.1 


4798 539 


Till 


1.08 


-2.75 


45.4 


-2.70 


44.7 


6161 295 


Cal 


2.52 


-1.18 


71.6 


-1.08 


69.7 


5211.544 


Till 


2.59 


-1.59 


32.1 


-1.54 


33.2 


6163 754 


Cal 


2.52 






-1.25 


83.4 


5336.783 


Till 


1.58 


-1.77 


71.5 


-1.63 


73.4 


6166 440 


Cal 


2.52 


-1.18 


71.5 


-1.02 


76.9 


5381.020 


Till 


1.57 


-1.91 


65.1 


-1.95 


59.5 


6169 044 


Cal 


2.52 


-0.75 


99.4 


-0.70 


97.0 


5418 756 


Till 


1.58 


-2.21 


49.7 


-2.17 


49.4 


6169 564 


Cal 


2.52 


-0.51 


117.6 


-0.53 


119.8 


5657.436 


VI 


1.06 








9.4 


6449 820 


Cal 


2.52 






-0.32 


127.8 


5668.362 


VI 


1.08 




8.7 




6.7 


6455 605 


Cal 


2.52 






-1.43 


53.7 


5670 851 


VI 


1.08 




21.1 




21.7 


6471 688 


Cal 


2.52 






-0.64 


101.5 


5727.661 


VI 


1.05 




10.9 




12.5 


6499 654 


Cal 


2.52 






-0.86 


87.2 


6090.216 


VI 


1.08 




35.6 




34.4 


4743.817 


Scl 


1.45 








8.5 


6135.370 


VI 


1.05 




12.6 




11.5 


5356 091 


Scl 


1.86 








1.8 


6150 154 


VI 


0.30 




12.6 




11.0 


5392 075 

^ .y y ^ t\j 1 ^ 


Scl 


1.99 








7.0 


6199 186 


VI 


0.29 




15.1 




14.0 


5484.611 


Scl 


1.85 








3.2 


6216.358 


VI 


0.28 








37.0 


5671 826 


Scl 


1.45 








19.3 


6274.658 


VI 


0.27 








8.7 


6239.408 


Scl 


0.00 








9.0 


6285.165 


VI 


0.28 




10.4 




16.2 


5318 346 


Sell 


1.36 








18.2 


4575 092 


CrI 


3.37 






-0.88 


13.6 


5357.190 


Sell 


1.51 




5.4 




5.2 


4616.120 


CrI 


0.98 






-1.31 


91.9 


5526.815 


Sell 


1.77 




77.9 




77.6 


4626.174 


CrI 


0.97 






-1.47 


84.8 


5657.874 


Sell 


1.51 




69.8 




69.2 


4708.019 


CrI 


3.17 






-0.06 


58.1 


5684.189 


Sell 


1.51 




40.8 




41.2 


4737.355 


CrI 


3.09 






-0.06 


62.0 


6245.660 


Sell 


1.51 




38.0 




35.6 


4756.137 


CrI 


3.10 


0.09 


74.6 


0.03 


66.4 


6320.867 


Sell 


1.50 




9.1 




8.4 


4801.047 


CrI 


3.12 


-0.28 


51.5 


-0.28 


46.2 


4518.023 


Til 


0.83 


-0.49 


76.2 






4936.335 


CrI 


3.11 


-0.35 


48.5 


-0.32 


47.2 


4548.765 


Til 


0.83 


-0.55 


73.7 






4964.916 


CrI 


0.94 






-2.50 


41.6 


4562.625 


Til 


0.02 






-2.73 


11.7 


5200.207 


CrI 


3.38 


-0.58 


24.6 


-0.50 


26.8 


4617.254 


Til 


1.75 






0.23 


64.6 


5214.144 


CrI 


3.37 


-0.77 


17.7 


-0.73 


18.4 


4758.120 


Til 


2.25 


0.44 


56.0 


0.26 


43.3 


5238.964 


CrI 


2.71 


-1.43 


16.6 


-1.36 


17.9 


4759.272 


Til 


2.25 


0.47 


58.0 


0.25 


45.9 


5247.566 


CrI 


0.96 


-1.73 


83.1 


-1.61 


85.1 


4778.259 


Til 


2.24 


-0.38 


18.0 


-0.38 


17.2 


5272.007 


CrI 


3.45 






-0.36 


30.3 
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Table 8. continued. 



A [A] 


Id. 


X 

[eV] 


Gany 1 


Gany 2 


A [A] 


Id. 


X 

[eV] 


Gany 1 


Gany 2 




logg/ 


EW 


logg/ 


EW 




logg/ 


EW 


logg/ 


EW 


5287.183 


CrI 


3.44 


-0.90 


12.0 


-0.86 


12.6 


5225.525 


Fel 


0.11 


-4.81 


11.1 


-4.72 


75.1 


5296.691 


CrI 


0.98 


— 


— 


-1.40 


93.5 


5242.491 


Fel 


3.63 


-1.26 


86.0 


-1.16 


Sl.l 


5300.751 


CrI 


0.98 


-2.13 


63.3 


-2.11 


56.7 


5243.773 


Fel 


4.26 


-1.13 


62.1 


-1.05 


63.8 


5304.183 


CrI 


3.46 


— 


— 


-0.72 


16.0 


5247.049 


Fel 


0.09 


-5.00 


70.4 


-4.93 


70.0 


5318.810 


CrI 


3.44 


-0.66 


19.3 


-0.65 


18.8 


5250.216 


Fel 


0.12 


-4.86 


75.4 


-4.95 


62.4 


5628.621 


CrI 


3.42 


— 


— 


-0.82 


14.6 


5320.040 


Fel 


3.64 


-2.47 


23.6 


-2.44 


23.5 


5648.279 


CrI 


3.82 


0.90 


5.7 


— 


— 


5321.109 


Fel 


4.43 


-1.24 


47.1 


-1.21 


45.8 


5784.976 


CrI 


3.32 


-0.45 


34.0 


-0.39 


33.8 


5332.908 


Fel 


1.56 


-3.07 


92.2 


-2.84 


95.7 


5787.965 


CrI 


3.32 


— 


— 


-0.12 


50.4 


5379.574 


Fel 


3.69 


-1.60 


64.9 


-1.56 


62.6 


6330.097 


CrI 


0.94 


-2.88 


30.0 


-2.90 


26.9 


5389.486 


Fel 


4.41 


-0.63 


84.4 


-0.56 


83.3 


4588.203 


CrII 


4.07 


— 


— 


-0.73 


71.7 


5395.222 


Fel 


4.44 


-1.73 


22.9 


-1.73 


21.8 


4592.049 


CrII 


4.07 


-1.30 


50.1 


-1.23 


50.8 


5412.791 


Fel 


4.43 


-1.76 


21.9 


-1.75 


21.4 


5305.855 


CrII 


3.83 


— 


— 


-2.06 


27.7 


5432.946 


Fel 


4.44 


-0.79 


72.7 


-0.69 


74.4 


5308.377 


CrII 


4.07 


— 


— 


-1.81 


27.1 


5436.297 


Fel 


4.39 


-1.36 


42.5 


-1.31 


42.8 


5313.526 


CrII 


4.07 


— 


— 


-1.61 


34.9 


5473.168 


Fel 


4.19 


-1.99 


21.9 


-1.96 


22.2 


5502.025 


CrII 


4.17 


— 


— 


-1.87 


21.8 


5483.108 


Fel 


4.15 


-1.46 


49.1 


-1.45 


47.1 


4626.538 


Mnl 


4.71 


— 


30.5 


— 


26.8 


5491.845 


Fel 


4.19 


-2.19 


15.2 


-2.23 


13.8 


4739.113 


Mnl 


2.94 


— 


61.6 


— 


63.3 


5494.474 


Fel 


4.07 


-1.89 


31.1 


-1.94 


27.6 


5004.892 


Mnl 


2.92 


— 


— 


— 


18.2 


5522.454 


Fel 


4.21 


-1.50 


44.2 


-1.44 


44.8 


5394.670 


Mnl 


0.00 


— 


— 


— 


81.8 


5560.207 


Fel 


4.43 


-1.12 


54.3 


-1.09 


53.0 


5399.479 


Mnl 


3.85 


— 


— 


— 


40.0 


5577.013 


Fel 


5.03 


-1.52 


11.9 


-1.49 


12.5 


5413.684 


Mnl 


3.86 


— 


24.9 


— 


25.7 


5587.573 


Fel 


4.14 


-1.54 


45.7 


-1.56 


41.8 


5420.350 


Mnl 


2.14 


— 


86.6 


— 


90.3 


5635.824 


Fel 


4.26 


-1.58 


37.8 


-1.55 


37.1 


5432.548 


Mnl 


0.00 


— 


52.3 


— 


53.8 


5636.705 


Fel 


3.64 


-2.51 


22.4 


-2.52 


21.1 


5537.765 


Mnl 


2.19 


— 


41.1 


— 


37.3 


5638.262 


Fel 


4.22 


-0.89 


78.8 


-0.79 


79.4 


6013.497 


Mnl 


3.07 


— 


87.0 


— 


86.9 


5641.436 


Fel 


4.26 


-1.04 


67.6 


-0.99 


66.1 


6021.803 


Mnl 


3.07 


— 


97.0 


— 


99.0 


5646.697 


Fel 


4.26 


-2.38 


9.2 


-2.48 


5.5 


4523.400 


Fel 


3.65 


-1.85 


51.4 


— 


— 


5650.019 


Fel 


5.10 


-0.78 


39.5 


-0.82 


35.8 


4537.676 


Fel 


3.27 


-2.96 


17.6 


— 


— 


5652.319 


Fel 


4.26 


-1.79 


27.6 


-1.79 


26.1 


4556.925 


Fel 


3.25 


-2.69 


28.8 


— 


— 


5661.348 


Fel 


4.28 


-1.88 


23.3 


-1.81 


24.6 


4585.343 


Fel 


4.61 


-1.57 


22.6 


-1.59 


21.1 


5680.240 


Fel 


4.19 


-2.34 


11.5 


-2.30 


12.1 


4593.555 


Fel 


3.94 


-2.03 


29.2 


-2.00 


29.1 


5701.557 


Fel 


2.56 


-2.20 


90.2 


-2.13 


86.6 


4598.125 


Fel 


3.28 


-1.61 


82.2 


-1.61 


76.3 


5705.473 


Fel 


4.30 


-1.35 


48.0 


-1.44 


40.8 


4602.000 


Fel 


1.61 


-3.32 


74.7 


-3.21 


73.1 


5731.761 


Fel 


4.26 


-1.13 


62.4 


-1.14 


58.0 


4741.535 


Fel 


2.83 


-2.21 


73.2 


-2.08 


73.8 


5738.240 


Fel 


4.22 


-2.12 


16.6 


-2.19 


14.2 


4749.961 


Fel 


4.56 


-1.17 


44.4 


-1.28 


36.2 


5775.069 


Fel 


4.22 


-1.21 


60.0 


-1.11 


62.1 


4793.961 


Fel 


3.05 


-3.40 


11.4 


-3.54 


8.5 


5778.463 


Fel 


2.59 


-3.53 


22.5 


-3.47 


23.7 


4794.355 


Fel 


2.42 


-3.88 


14.8 


-3.86 


14.5 


5784.666 


Fel 


3.40 


-2.59 


29.1 


-2.53 


29.4 


4798.273 


Fel 


4.19 


-1.48 


45.4 


-1.43 


45.5 


5811.916 


Fel 


4.14 


-2.40 


11.2 


-2.36 


11.8 


4798.743 


Fel 


1.61 


-4.19 


33.5 


-4.22 


35.1 


5814.805 


Fel 


4.28 


-1.85 


24.6 


-1.85 


23.4 


4808.147 


Fel 


3.25 


-2.65 


31.2 


-2.59 


33.8 


5835.098 


Fel 


4.26 


-2.11 


16.0 


-2.10 


15.7 


4907.733 


Fel 


3.43 


-1.80 


65.8 


-1.76 


63.5 


5849.681 


Fel 


3.69 


-2.98 


8.3 


-2.90 


9.5 


4908.032 


Fel 


3.93 


-1.56 


39.7 


-1.77 


40.2 


5852.222 


Fel 


4.55 


-1.26 


40.5 


-1.16 


43.2 


4911.788 


Fel 


3.93 


-1.72 


45.6 


-1.61 


48.3 


5855.086 


Fel 


4.61 


-1.58 


22.9 


-1.49 


25.3 


4961.915 


Fel 


3.63 


-2.38 


26.9 


-2.31 


28.6 


5856.096 


Fel 


4.29 


-1.64 


33.7 


-1.52 


37.3 


4962.565 


Fel 


4.18 


-1.33 


54.5 


-1.25 


55.3 


5859.596 


Fel 


4.55 


-0.70 


72.7 


-0.60 


74.4 


4969.916 


Fel 


4.22 


-0.89 


77.6 


-0.78 


79.5 


5916.249 


Fel 


2.45 


-2.97 


57.0 


-2.89 


56.8 


5023.189 


Fel 


4.28 


-1.47 


41.7 


-1.37 


44.3 


5927.786 


Fel 


4.65 


-1.12 


43.0 


-1.05 


44.5 


5025.091 


Fel 


4.26 


-1.87 


23.6 


-1.83 


24.1 


5929.666 


Fel 


4.55 


-1.17 


45.4 


-1.18 


42.4 


5025.313 


Fel 


4.28 


-1.97 


19.2 


-1.78 


25.3 


5930.173 


Fel 


4.65 


-0.33 


91.9 


-0.26 


91.2 


5054.647 


Fel 


3.64 


o no 


4U.0 


1 m 

— i.y/ 


AAA 

44.4 


5956.692 


Fel 


0.86 


A 

—4.03 


A 


— 4.j4 




5067.162 


Fel 


4.22 


-0.98 


72.9 


-0.90 


74.7 


6005.551 


Fel 


2.59 


-3.50 


24.2 


-3.48 


23.6 


jUIZ.OI 1 


Fel 


4.22 


-1.09 


66.0 


-0.98 


67.9 


6007.968 


Fel 


4.65 


-0.73 


66.0 


-0.71 


63.3 


5109.649 


Fel 


4.30 


-0.77 


81.1 


-0.68 


83.3 


6012.212 


Fel 


2.22 


-3.82 


26.0 


-3.79 


25.5 


5127.359 


Fel 


0.93 


-3.57 


96.3 


-3.30 


101.1 


6078.499 


Fel 


4.79 


-0.36 


82.1 


-0.29 


81.9 


5127.680 


Fel 


0.05 


-5.97 


24.3 


-5.84 


27.8 


6079.014 


Fel 


4.65 


-1.04 


47.6 


-0.98 


48.6 


5196.065 


Fel 


4.26 


-0.90 


75.4 


-0.78 


77.3 


6082.708 


Fel 


2.22 


-3.58 


37.8 


-3.53 


37.3 


5197.929 


Fel 


4.30 


-1.52 


38.5 


-1.50 


37.0 


6093.666 


Fel 


4.61 


-1.37 


32.4 


-1.34 


31.6 


5213.818 


Fel 


3.94 


-2.75 


7.7 


-2.67 


8.9 


6098.250 


Fel 


4.56 






-1.75 


19.5 


5223.188 


Fel 


3.63 


-2.29 


31.4 


-2.26 


31.0 


6120.249 


Fel 


0.92 


-5.81 


6.7 


-5.86 


5.8 
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Table 8. continued. 



A [A] 


Id. 


X 

[eV] 


Gany 1 


Gany 2 


A [A] 


Id. 


X 

[eV] 


Gany 1 


Gany 2 




logs/ 


EW 


logg/ 


EW 




logg/ 


EW 


logg/ 


EW 


6137.002 


Fel 


2.20 


-2.93 


72.4 


-2.83 


71.2 


6455.001 


Col 


3.63 


— 


— 


— 


14.4 


6151.616 


Pel 


2.18 


-3.35 


51.7 


-3.30 


50.6 


4935.831 


Nil 


3.94 


-0.41 


64.4 


-0.37 


72.0 


6159.382 


Fel 


4.61 


-1.86 


13.7 


-1.84 


13.7 


4946.029 


Nil 


3.80 


-1.18 


30.5 


-1.22 


27.5 


6173.340 


Fel 


2.22 


-2.95 


70.3 


-2.93 


66.1 


4953.200 


Nil 


3.74 


-0.82 


51.8 


-0.67 


56.1 


6187.987 


Fel 


3.94 


-1.71 


48.1 


-1.64 


48.4 


5010.934 


Nil 


3.63 


-0.92 


51.8 


-0.91 


49.1 


6199.508 


Fel 


2.56 


-4.25 


5.8 


-4.34 


4.3 


5032.723 


Nil 


3.90 


-1.20 


25.6 


-1.09 


28.8 


6200.321 


Fel 


2.61 


-2.45 


76.4 


-2.37 


74.6 


5094.406 


Nil 


3.83 


-1.14 


31.3 


-1.06 


33.1 


6213.428 


Fel 


2.22 


-2.68 


84.4 


-2.59 


82.1 


5197.157 


Nil 


3.90 


-1.15 


27.8 


-1.09 


28.8 


6219.287 


Fel 


2.20 


-2.52 


93.3 


-2.47 


88.3 


5220.300 


Nil 


3.74 


-1.30 


27.6 


-1.23 


29.4 


6226.730 


Fel 


3.88 


-2.10 


31.0 


-2.08 


29.5 


5392.330 


Nil 


4.15 


-1.24 


15.7 


-1.31 


13.3 


6240.645 


Fel 


2.22 


-3.37 


49.2 


-3.34 


46.6 


5435.866 


Nil 


1.99 


-2.47 


54.2 


-2.38 


54.1 


6265.131 


Fel 


2.18 


-2.65 


87.9 


-2.53 


86.8 


5452.860 


Nil 


3.84 


-1.48 


17.5 


-1.48 


18.5 


6271.283 


Fel 


3.33 


-2.70 


27.8 


-2.67 


27.8 


5494.876 


Nil 


4.10 


-1.11 


21.7 


-1.07 


22.7 


6297.792 


Fel 


2.22 


-2.77 


79.9 


-2.34 


96.5 


5587.853 


Nil 


1.93 


-2.37 


62.4 


-2.32 


60.4 


6315.813 


Fel 


4.07 


-1.66 


43.9 


-1.67 


41.2 


5625.312 


Nil 


4.09 


-0.59 


47.6 


-0.63 


42.9 


6322.691 


Fel 


2.59 


-2.47 


77.0 


-2.31 


79.0 


5628.354 


Nil 


4.09 


— 


— 


-1.28 


15.8 


6358.687 


Fel 


0.86 


— 


— 


-3.79 


87.8 


5637.128 


Nil 


4.09 


-0.80 


36.1 


-0.80 


34.8 


6380.750 


Fel 


4.19 


— 


— 


-1.30 


53.9 


5748.346 


Nil 


1.68 


-3.26 


30.7 


-3.22 


30.5 


6385.726 


Fel 


4.73 


— 


— 


-1.94 


9.1 


5846.986 


Nil 


1.68 


-3.40 


24.7 


-3.33 


26.3 


6392.538 


Fel 


2.28 


— 


— 


-3.98 


17.3 


6086.276 


Nil 


4.26 


-0.44 


47.7 


-0.44 


45.4 


6393.612 


Fel 


2.43 


— 


— 


-1.60 


132.8 


6176.807 


Nil 


4.09 


-0.26 


66.9 


-0.28 


61.6 


6430.856 


Fel 


2.18 


— 


— 


-2.01 


117.8 


6177.236 


Nil 


1.83 


-3.44 


18.4 


-3.52 


14.8 


6498.945 


Fel 


0.96 


— 


— 


-4.58 


48.4 


6186.709 


Nil 


4.10 


-0.90 


31.6 


-0.87 


31.6 


4576.339 


Fell 


2.84 


-3.13 


66.0 


-3.03 


66.3 


6191.187 


Nil 


1.68 


-2.30 


80.2 


-2.18 


79.3 


4656.981 


Fell 


2.89 


— 


— 


-3.59 


41.2 


6327.604 


Nil 


1.68 


-3.07 


41.2 


-3.04 


39.8 


4720.149 


Fell 


3.20 


-4.49 


6.0 


-4.57 


6.0 


6370.357 


Nil 


3.54 


— 


— 


-1.75 


18.5 


4993.358 


Fell 


2.81 


-3.74 


40.0 


-3.69 


40.8 


6378.256 


Nil 


4.15 


— 


— 


-0.77 


34.3 


5197.576 


Fell 


3.23 


-2.45 


81.4 


-2.32 


82.3 


5218.209 


Cul 


3.82 


— 


54.1 


— 


55.8 


5234.630 


Fell 


3.22 


-2.36 


86.1 


-2.23 


86.8 


5220.086 


Cul 


3.82 


— 


17.3 


— 


17.2 


5264.812 


Fell 


3.33 


-2.98 


52.5 


-2.96 


50.7 


5782.136 


Cul 


1.64 


— 


79.9 




85.0 


5325.560 


Fell 


3.22 


-3.16 


49.0 


-3.15 


47.3 


4810.537 


Znl 


4.08 


-0.33 


76.1 


-0.27 


79.0 


5414.075 


Fell 


3.22 


-3.60 


29.1 


-3.54 


31.0 


4607.338 


Sri 


0.00 


0.02 


48.8 


0.12 


48.1 


5425.257 


Fell 


3.20 


-3.23 


46.5 


-3.25 


43.6 


4883.690 


YII 


1.08 


— 


— 


0.06 


62.8 


5427.826 


Fell 


6.72 


— 




-1.31 


6.7 


4900.124 


YII 


1.03 


-0.29 


55.0 


-0.07 


60.1 


6084. 1 1 1 


Fell 


3.20 


-3.75 


24. 1 


-3.78 


22.7 


5087.426 


YII 


1.08 


-0.43 


46.9 


-0.33 


48.5 


6149.249 


Fell 


3.89 


-2.76 


38.5 


-2.73 


38.6 


5200.415 


YII 


0.99 


-0.70 


39.0 


-0.71 


36.5 


6247.562 


Fell 


3.89 


-2.37 


56.4 


-2.37 


53.9 


5289.820 


YII 


1.03 


-1.77 


5.4 


— 


— 


6369.463 


Fell 


2.89 


-4.14 


20.9 


-4.15 


20.3 


5402.780 


YII 


1.84 


-0.48 


15.6 


-0.61 


15.0 


6383.715 


Fell 


5.55 


— 


— 


-2.07 


10.7 


4739.454 


Zrl 


0.65 


0.00 


7.5 


0.04 


7.9 


6385.458 


Fell 


5.55 


— 


— 


-2.44 


5.1 


4613.921 


Zrll 


0.97 


-0.62 


37.1 


-0.61 


35.7 


6416.928 


Fell 


3.89 


-2.65 


43.5 


-2.69 


40.5 


5112.279 


Zrll 


1.66 


-0.75 


10.2 


-0.81 


12.3 


6456.391 


Fell 


3.90 


— 


— 


-2.24 


59.4 


5853.688 


Ball 


0.60 


-1.01 


64.0 


-0.84 


65.6 


4749.662 


Col 


3.05 


— 


50.0 


— 


40.0 


6141.727 


Ball 


0.70 


0.13 


120.0 


0.24 


119.6 


4792.862 


Col 


3.25 










6496.908 


Ball 


0.60 


r\ nc 
— U.Uj 


lUl.U 


— U.U / 


lUD.Z 


4813.479 


Col 


3.21 




48.0 




48.6 


4523.080 


Cell 


0.52 


0.24 


21.4 






JZ1/.091 


Co i 


TCI 
3.51 




18.4 




20.4 


450/. 36/ 


TT 

Ceil 


0.48 


0.37 


27.9 






5280.629 


Col 


3.63 








20.6 


4628.160 


Cell 


0.52 


0.21 


20.4 


0.27 


21.8 


5342.708 


Col 


4.02 




32.1 




32.1 


4773.959 


Cell 


0.92 


0.32 


12.7 


0.31 


12.1 


5359.192 


Col 


4.15 




11.0 




10.4 


5274.236 


Cell 


1.04 


0.48 


14.3 


0.40 


12.0 


5381.772 


Col 


4.24 




9.9 




6.0 


5089.831 


Ndll 


0.20 


-1.31 


2.7 


-1.23 


3.5 


5454.572 


Col 


4.07 








18.2 


5319.820 


Ndll 


0.55 


-0.21 


14.6 


-0.17 


15.1 


5647.234 


Col 


2.28 




15.1 




14.8 


4566.233 


Smll 


0.33 


-0.19 


12.3 


-0.19 


12.0 


6000.678 


Col 


3.62 




4.3 




5.8 
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Table 9. Oscillator strengths (log gf) for lines with hyperfine structure computed based on the Ganymede spectra observed in the 
first (Gany 1) and second (Gany 2) runs. 





logs/ 
Gany 1 Gany 2 


) lAl 


log gf 
Gany 1 Gany 2 


) lAl 


Gany 1 Gany 2 


1 lAl 


logs/ 
Gany 1 Gany 2 


) lAl 


logs/ 
Gany 1 Gany 2 


Mgl 


: 4571.102 




Sc 11 


: 5657.874 




V 1 : 


6274.658 




Mn 1 


: 5537.765 




Co 1 


: 5381.772 




4j / l.U/o 


-6.67 


-6.49 


565 /.808 


-1.33 


-1.32 


£OTy| <;^A 

6z/4.64U 




-2.09 


5537.691 


-2.76 


-2.79 


CO O 1 AAC 

5381.695 


-0.47 


-0.69 


4 J / 1 .Uo / 


-6.71 
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-1.22 


ATT/1 ACQ 
Oz /4.0j8 




-2.09 


CCOO 010 
JJJ / . / lU 


-2.64 
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CO 1 OOO 
J JSI . / J8 


-0.55 


-0.77 


4 J / 1 .uyo 


-5.81 


-5.63 
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-1.44 
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Oz /4.0 /O 




-2.09 


CCOO ony 

JJJ /. /y8 


-2.61 


-2.64 


CO 1 OOO 
J J81 . / /Z 


-0.63 


-0.85 


Mgl 


: 4730.038 




jOj /.oyj 


-1.33 


-1.32 


V 1 : 


6285.165 




CCOO OA/1 
JJJ /. /04 


-2.64 


-2.67 


COQl OOO 

J j8i . /yy 


-0.74 


-0.96 


4/jU.Ujl 


-3.13 


-3.10 


Sc 11 


: 5684.189 




ATQ C } An 
0Z8j. 14 / 


-2.11 


-1.86 


CCOO QAO 
JJJ /.8Uz 


-2.28 


-2.31 


COQl QO/1 
J J81 .8Z4 
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^TJ A Al O 

4/ jU.Ujo 


-3.17 


-3.14 


5684.123 


-1.72 


-1.71 


6285.165 


-2.11 


-1.86 


Mnl 


: 6013.497 




Col 


: 5454.572 




4 / JU.U40 


-2.27 


-2.24 


5684. 156 


-1.62 


-1.61 


ATQC 1 QO 
OzSj. 18j 


-2.11 


-1.86 


AO 10 /I 0/1 
OUl J.4 /4 


-0.75 


-0.67 


C/1C/1 /IOC 
J4j4.4yj 




-0.38 


Mgl 


: 5711.095 




CAQ/1 1 OA 

j064. lay 


-1.83 


-1.82 


Mnl 


: 4626.538 




AO 1 /I QA 
OUl J. 480 


-0.96 


-0.88 


C/1 C/1 COQ 
j4j4.jj8 




-0.45 


J / i 1 .U /4 


-2.69 


-2.69 


c/;c/i OAC 
j0o4.ZU6 


-1.72 


-1.71 


/I AOA /I A/1 
40Z0.404 


-0.94 


-0.99 


AO 1 COl 
OUl J.jUl 


-1.10 


-1.02 


5454.572 




-0.54 


J / i 1 .US J 


-2.73 


-2.73 


Sc 11 


: 6245.660 




/1A0A CA/I 
40Z0.JU4 


-0.14 


-0.19 


AO 1 C 1 O 

OUl J.J ly 


-0.77 


-0.69 


C/1C/1 coo 
j4j4.jyy 




-0.64 


C"7 1 1 An 1 

J / 1 i.uyi 


-1.83 


-1.83 


6245.661 


-1.78 


-1.81 


4626.530 


-0.39 


-0.44 


AA 1 O COO 

6013.53/ 


-1.35 


-1.27 


C /I C /I AO /I 

5454.624 




-0.59 


Sc I 


4743.817 




0z4j.04z 


-1.89 


-1.92 


4626.565 


-0.19 


-0.24 


Mn 1 


: 6021.803 




Co 1 


: 5647.234 




4743.75 1 




-0.37 


AO/IC *^AA 

Oz4j.ouy 


-1.68 


-1.71 


/1ATA CTO 
40Z0.J / J 


-0.49 


-0.54 


AOO 1 OA/1 

OUZl . /04 


-1.37 


-1.24 


5647. 191 


-2.18 


-2.17 


4 /4J. /54 




-0.28 


0Z4J.J /O 


-1.77 


-1.80 


Mnl 


: 4739.113 




AOOl OCO 
OUZ 1 . / 8U 


-1.22 


-1.09 


eA/lO OAO 
J04/.ZUU 


-2.25 


-2.24 


4 /4J.O 1 / 




-0.49 


Sc 11 


: 6320.867 




4 / jy.uyy 


-1.30 


-1.23 


AOOl OOO 

OUZl . /y / 


-0.39 


-0.26 


CA/IO 00/1 
j04/.Zj4 


-2.34 


-2.33 


AnAl 

4 /4J.0JO 




-0.38 


OjZU.554 


-2.52 


-2.57 


/1"70O 110 

4/ jy. 1 1 J 


-1.44 


-1.37 


AOO 1 QAA 
OUZl .8uO 


-0.60 


-0.47 


5647.261 


-2.44 


-2.43 


Sc I 


5356.091 




*^00A Q^iC 
OjZU.OOJ 


-2.63 


-2.68 


/l"70O 1 OA 

4 / jy. izo 


-1.60 


-1.53 


AOO 1 Q 1 /I 
OUZl .8 14 


-0.47 


-0.34 


CA/IO OO 1 

j04/.zyi 


-2.39 


-2.38 






-0.67 


OjZU.ojz 


-2.42 


-2.47 


/1"70O 1 /I C 

4/ jy. i4j 


-1.15 


-1.08 


Co 1 


: 4749.662 




Co 1 


: 6000.678 




J JJO.UJO 




-0.58 


iiOOA TAA 

ojzu. /yy 


-2.52 


-2.57 


/l"70O 1 A"7 

4 / jy. 10 / 


-2.50 


-2.43 


/lO/IO A/1 1 
4/4y.041 


-0.63 


-0.77 


AOAO AAO 
OUUU.OU / 


-2.21 


-2.09 


jjjo.uyi 




-0.79 


V 1 


5657.436 




Mn 1 


: 5004.892 




/lO/lO AOC 
4 /4y.O / J 


-0.81 


-0.85 


AOAO A/IO 

ouuu.o4y 


-1.50 


-1.38 


5356.110 




-0.68 


5657.418 




-1.34 


CAA^ OTO 

5U04.8 /8 




-2.08 


^O ^ A OA/1 

4/4y. /04 


-0.92 


-1.06 


A AAA AOO 
6U0U.6/8 


-1.38 


-1.26 


Sc I 


5392.075 




5657.436 




-1.34 


cnci/i om 
jUU4.oyz 




-2.22 


/lO/lO OOO 

4 /4y. /zy 


-0.88 


-1.02 


AOAO O 1 O 
OUUU. / 1 / 


-1.55 


-1.43 


com nno 




0.04 


5657.454 




-1.34 


CAA/1 one 
jUU4.yUj 




-2.38 


/lO/lO OA/1 
4/4y. /04 


-2.28 


-2.42 


AOAO OCO 
OUUU. / JZ 


-1.54 


-1.42 


com n^n 




0.13 


V 1 


5668.362 




cnci/i O0/1 
jUU4.yz4 




-1.93 


Co 1 


: 4792.862 




Co 1 


: 6455.001 




com mc 




-0.08 


5668.344 


-1.37 


-1.47 


CAA/1 A^A 

5U04.y46 




-3.28 


/I OAO O 1 1 

4/yz.8ii 


-2.12 


-2.11 


A /I c /I AO 1 

6454. y31 




-0.56 


CO no no .1 
53yz.U64 




0.03 


5668.362 


-1.37 


-1.47 


Mnl 


: 5394.670 




/lOAO OOO 

4/yZ.8Z/ 


-1.46 


-1.45 


A /I C /I AOA 

6454. y/y 




-1.35 


Sc I 


5484.611 




CA^iQ OOA 
J0o8. JSU 


-1.37 


-1.47 


COO/1 An 
jjy4.oi / 




-3.93 
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5671.845 
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Col 
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-1.06 


6216.376 




-1.38 


5432.598 


-4.80 


-4.74 


5359.115 


-0.51 


-0.53 


5782.113 


-2.84 


-2.74 


5526.834 


-0.97 


-0.95 














5359.158 


-0.59 


-0.61 


5782.124 


-2.84 


-2.74 




















5359.192 


-0.67 


-0.69 


5782.153 


-2.74 


-2.64 




















5359.219 


-0.78 


-0.80 


5782.173 


-2.39 


-2.29 




















5359.244 


-0.73 


-0.75 
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Table 10. Elemental abundance relative to iron. Our results for carbon abundances are presented in Table |6] 



Star 


[Na/Fe] 


iMg/Fel 


[Si/Fe] 


[Ca/Fe] 


[Sc/Fe] 


[Ti/Fe] 


[V/Fe] 


[Cr/Fe] 


[Mn/Fe] 


[Co/Fe] 


HD 1835 


0.04 


0.00 


-0.01 


0.08 


-0.01 


-0.02 


0.08 


0.04 


0.04 


-0.03 


HD 20807 


0.05 


0.02 


0.00 


0.00 


0.06 


0.04 


0.04 


-0.02 


-0.06 


0.05 


HD 26491 


-0.02 


0.06 


0.02 


0.03 


0.06 


0.04 


0.07 


0.02 


0.00 


0.03 


HD 33021 


-0.05 


0.14 


0.08 


0.12 


0.11 


0.16 


0.11 


-0.01 


-0.18 


0.06 


HD 39587 


-0.10 


-0.01 


-0.04 


0.09 


0.01 


0.03 


-0.03 


0.03 


-0.06 


-0.05 


HD 43834 


0.00 


0.05 


0.00 


-0.02 


0.02 


-0.02 


0.04 


-0.01 


0.07 


0.02 


HD 50806 


0.04 


0.21 


0.09 


0.05 


0.06 


0.09 


0.08 


0.00 


-0.04 


0.08 


HD 53705 


0.03 


0.16 


0.06 


0.05 


0.13 


0.11 


0.07 


0.02 


-0.12 


0.09 


HD 84117 


0.08 


-0.03 


0.03 


-0.03 


-0.01 


-0.05 


0.03 


-0.02 


-0.09 


-0.04 


HD 102365 


- 


0.09 


0.03 


0.09 


0.06 


0.11 


0.05 


-0.02 


-0.11 


0.05 


HD 112164 


0.08 


0.07 


0.04 


-0.07 


-0.04 


-0.04 


0.05 


-0.01 


0.06 


0.05 


HD 114613 


0.08 


0.04 


0.02 


-0.02 


-0.01 


0.00 


0.02 


0.00 


0.08 


0.01 


HD 115383 


-0.08 


0.01 


-0.03 


-0.03 


0.02 


0.00 


0.04 


-0.03 


-0.01 


-0.02 


HD 115617 


-0.04 


0.06 


0.00 


0.00 


0.05 


0.02 


0.01 


-0.01 


0.03 


-0.02 


HD 117176 


-0.10 


0.07 


0.01 


0.02 


0.03 


0.01 


0.04 


-0.01 


-0.05 


0.02 


HD 128620 


0.06 


0.04 


0.01 


-0.06 


0.03 


0.01 


0.07 


-0.01 


0.07 


0.07 


HD 141004 


-0.05 


0.02 


-0.01 


-0.01 


0.04 


0.03 


0.04 


-0.03 


-0.04 


0.06 


1 1 1 J i40Zjj 


— U.Uj 


U.UU 


n (YX 

U.Uj 


n no 
u.uz 


u.u / 


— u.ut 


U.U/ 


U.Ui 


—U.Ui 


— U.U2 


HD 147513 




-0.06 


-0.04 


0.03 


0.05 


0.07 


0.12 


-0.02 


-0.07 


-0.03 


HD 160691 


0.09 


0.07 


0.03 


-0.06 


0.00 


-0.02 


0.07 


-0.03 


0.07 


0.03 


HD 177565 


0.02 


0.04 


0.01 


-0.01 


-0.04 


0.01 


0.03 


0.02 


0.05 


0.01 


HD 181321 


— 


-0.06 


0.03 


0.03 


-0.07 


0.01 


— 


0.03 


-0.13 


0.11 


HD 188376 


-0.03 


0.03 


0.02 


0.03 


-0.07 


-0.07 


0.03 


0.00 


0.01 


-0.06 


HD 189567 


-0.03 


0.11 


0.04 


0.00 


0.10 


0.11 


0.03 


0.01 


-0.12 


0.03 


HD 196761 


0.00 


0.06 


0.02 


-0.01 


0.06 


0.05 


-0.01 


0.01 


-0.02 


0.05 




[Ni/Fe] 


[Cu/Fe] 


[Zn/Fe] 


[Sr/Fe] 


[Y/Fe] 


[Zr/Fe] 


[Ba/Fe] 


[Ce/Fe] 


[Nd/Fe] 


[Sm/Fe] 


HD 1835 


-0.02 


-0.01 


0.07 


0.15 


0.02 


0.09 


0.19 


0.02 


- 


-0.29 


HD 20807 


-0.03 


0.02 


-0.03 


0.04 


-0.08 


0.04 


-0.06 


0.07 


0.02 


0.15 


HD 26491 


0.00 


0.03 


0.02 


-0.08 


-0.08 


-0.07 


0.03 


-0.06 


-0.11 


-0.08 


HD 33021 


-0.02 


-0.03 


-0.03 


-0.06 


-0.13 


0.13 


0.02 


-0.03 


0.06 


-0.02 


HD 39587 


-0.08 


-0.11 


-0.01 


0.19 


0.06 


0.16 


0.33 


0.16 


- 


-0.14 


HD 43834 


0.03 


0.03 


0.02 


0.06 


-0.04 


0.07 


-0.03 


0.00 


-0.06 


0.06 


HD 50806 


0.04 


0.15 


0.10 


-0.08 


-0.06 


0.11 


0.00 


0.01 


-0.09 


0.10 


HD 53705 


-0.02 


0.05 


0.13 


-0.09 


-0.15 


0.00 


-0.02 


-0.08 


-0.07 


-0.13 


HD 84117 


-0.04 


-0.03 


-0.03 


0.09 


0.07 


0.15 


0.15 


-0.02 


0.09 


0.00 


HD 102365 


-0.02 


0.02 


0.07 


-0.02 


-0.06 


-0.10 


0.12 


-0.06 


0.04 


0.02 


HD 112164 


0.03 


0.14 


0.00 


-0.06 


-0.16 


-0.06 


-0.03 


-0.09 


-0.09 


-0.10 


HD 114613 


0.03 


0.10 


0.02 


-0.02 


-0.06 


-0.01 


0.00 


-0.14 


-0.11 


-0.02 


HD 115383 


-0.03 


0.05 


-0.07 


0.00 


0.01 


0.06 


0.20 


-0.04 


0.16 


0.02 


HD 115617 


-0.01 


-0.01 


0.08 


0.05 


-0.11 


0.01 


-0.02 


0.05 


-0.07 


0.28 


HD 117176 


0.00 


0.07 


-0.06 


-0.03 


-0.06 


0.04 


0.06 


0.06 


0.19 


0.11 


HD 128620 


0.04 


0.14 


-0.08 


- 


-0.06 


0.06 


-0.08 


-0.06 


-0.01 


0.04 


HD 141004 


-0.01 


0.10 


-0.05 


-0.02 


-0.02 


-0.05 


0.03 


-0.07 


- 


-0.10 


HD 146233 


-0.01 


-0.02 


-0.04 


0.17 


0.00 


0.06 


0.12 


0.02 


0.13 


0.09 


HD 147513 


-0.06 


-0.11 


-0.08 




0.22 


0.24 


0.43 


0.06 


0.18 




HD 160691 


0.03 


0.18 


0.11 


0.01 


-0.06 


0.07 


0.01 


-0.08 


-0.09 


0.09 


HD 177565 


0.04 


0.04 


0.02 


-0.06 


-0.04 


-0.09 


-0.01 


0.02 


-0.10 


0.07 


HD 181321 


-0.09 




-0.14 




-0.04 




0.23 








HD 188376 


-0.02 


-0.03 


-0.10 




0.00 


0.04 


0.21 


-0.07 


0.20 


-0.19 


HD 189567 


0.00 


-0.04 


0.03 


-0.07 


-0.07 


0.03 


0.06 


0.08 


0.06 


0.14 


HD 196761 


0.00 


0.01 


0.08 


-0.09 


-0.16 


0.05 


-0.03 


-0.04 


-0.04 


0.26 



